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Abstract 


This  work  is  a  theoretical  study  of  microwave  photonic  links  which  use  optical  fre¬ 
quency  modulation  (FM)  and  filter-slope  discrimination  for  demodulation.  The  high 
modulation  efficiency  of  optical  FM  devices  is  attractive  for  achieving  low  noise-figure 
finks,  but  linear  demodulation  of  the  signals  is  also  desired.  In  order  to  design  dis¬ 
criminators  which  produce  low  distortion,  this  paper  presents  general  full-signal  and 
small-signal  models  of  the  effects  of  arbitrary  optical  filtering  on  an  FM  fink,  includ¬ 
ing  the  interaction  between  FM  and  residual  IM.  The  small  signal  model  is  used  to 
derive  figures  of  merit  for  the  linearity,  noise  and  dynamic  range  of  the  FM  finks. 
The  results  of  the  models  invalidate  a  common  assumption:  that  linear  FM  to  IM 
discrimination  is  possible.  Instead,  the  discrimination  must  be  from  FM  to  amplitude 
modulation  of  the  electric  field  for  the  fink  response  to  be  linear. 

The  linearity  of  links  using  two  different  sets  of  FIR  discriminator  filters,  designed 
using  the  minimax  relative  error  and  maximally  linear  criteria,  are  compared  in  order 
to  evaluate  the  better  design  method.  The  evolution  of  fink  linearity  with  filter  order 
is  studied.  The  analytical  models  are  compared  to  numerical  simulations  of  the  filters, 
and  the  models  are  found  to  be  consistent.  Finally,  a  Monte  Carlo  simulation  is  used 
to  analyze  the  sensivity  to  errors  in  fabrication  for  high-order  filters  realized  in  planar 
lightwave  circuits  (PLC)  with  a  lattice  filter  architecture. 
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Chapter  1 


Introduction 

1.1  Microwave  photonic  links 

Microwave  photonics  is  the  application  of  photonic  devices,  such  as  lasers,  to  mi¬ 
crowave  systems  for  signal  generation,  signal  transmission  and  signal  processing.  Mi¬ 
crowave  photonic  systems  are  analog  in  the  sense  that  they  manipulate  arbitrary 
baseband  signals  as  well  as  digital  signals  that  are  modulated  onto  a  higher  carrier 
frequency.  Microwave  photonic  links  are  fiber  optic  links  which  transmit  microwave 
signals  to  remote  locations. 

Microwave  photonic  finks  have  been  explored  for  replacing  traditional  coaxial  links 
in  a  variety  of  applications  because  of  their  advantages  in  size,  weight,  immunity  to 
electromagnetic  interference,  bandwidth  and  power  consumption  [1,  2,  3].  The  most 
successful  commercial  applications  have  been  in  hybrid-fiber-coax  (HFC)  infrastruc¬ 
ture  for  distributing  cable-television  signals  and  in  hybrid-fiber-radio  (HFR)  for  dis¬ 
tributing  cellular  signals  to  remote  antennas  [1,  4].  Military  radar  and  communication 
systems  also  use  analog  fiber  optic  systems  for  antenna  remoting. 

The  applications  of  microwave  photonic  links  are  limited  by  the  noise  and  dis¬ 
tortion  performance  required  to  receive  the  transmitted  signals.  The  most  common 
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microwave  photonic  links  are  intensity-modulation  direct-detection  links  (IM-DD), 
where  the  optical  power  is  varied  in  proportion  to  the  input  signal,  transmitted  over 
optical  fiber  and  finally  directly-detected  by  a  photodetector.  The  noise  and  dis¬ 
tortion  performance  of  these  finks,  as  reviewed  by  [5],  are  not  yet  competitive  with 
electronic  systems  and  are  unsuitable  for  next  generation  military  wireless  systems. 
IM-DD  finks  have  high  (-30  dB  to  -20  dB)  RF-to-RF  signal  loss  due  to  low  modulation 
efficiency,  which  contributes  greatly  to  poor  noise  figures.  This  provides  an  impetus 
to  study  other  fink  architectures. 

1.2  Frequency  modulation 

Modulation  is  not  limited  to  the  intensity,  as  other  parameters  of  the  light  can  be 
used  to  convey  information.  The  amplitude,  phase,  frequency  and  polarization  of  the 
fight’s  electric  field  are  some  of  these  quantities  that  can  be  modulated.  Frequency 
modulated  (FM)  links,  where  the  instaneous  optical  frequency  of  the  laser  is  varied 
with  the  input  signal,  are  considered  to  be  promising  alternatives  to  IM  links.  Directly 
modulated  FM  lasers  have  been  demonstrated  with  high  modulation  efficiency  and 
with  modulation  bandwidths  that  are  not  limited  by  the  laser  relaxation  frequency  [6]. 
Recent  work  on  multi-section  DFB  [7]  and  EML  lasers  [8]  have  produced  modulation 
efficiencies  two  orders  of  magnitude  better  than  traditional  intensity  modulation. 
An  improvement  in  modulation  efficiency  could  make  an  major  impact  on  the  noise 
performance  of  microwave  photonic  finks. 

1.3  FM  discriminators 

Because  photodetectors  respond  to  the  intensity  envelope  of  the  fight  and  not  the 
instantaneous  frequency,  a  number  of  approaches  have  been  developed  to  demodulate 
the  FM  at  the  output  of  the  fink.  These  include  homodyne  detection,  heterodyne  de- 
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tection  and  discriminator  (or  slope)  detection.  The  first  two  systems  are  considered 
coherent  detection  systems,  since  two  beams  are  combined  coherently  at  the  pho¬ 
todetector.  In  the  third  approach,  an  optical  filter  acts  as  an  FM  discriminator  that 
converts  FM  to  amplitude  modulation  (AM),  and  the  AM  is  directly  detected  by  the 
photodetector.  Following  the  nomenclature  of  [9],  such  a  system  with  a  discriminator 
is  termed  an  FM  direct-detection  (FM-DD)  fink.  The  analysis  of  FM-DD  finks  is  the 
main  focus  of  this  work. 

FM-DD  links  use  optical  filters  as  discriminators  to  convert  FM  to  AM.  Paraphras¬ 
ing  [10],  the  sidebands  of  a  frequency-modulated  signal  possess  certain  amplitude  and 
phase  relationships  among  themselves  such  that  the  envelope  of  the  signal  is  inde¬ 
pendent  of  time.  A  discriminator  works  by  modifying  these  phase  and  amplitude 
relationships  such  that  the  amplitude  of  the  envelope  of  the  resultant  signal  fiuctu- 
ates  in  the  same  manner  versus  time  as  did  the  instantaneous  frequency  of  the  original 
signal. 

Heuristically,  one  can  think  of  the  FM  discriminator  as  a  function  with  a  frequency 
dependent  amplitude.  The  slope  of  the  function  converts  variations  in  the  optical 
frequency  into  variations  in  the  amplitude.  This  view  is  accurate  for  slow  variations 
of  the  optical  frequency.  However,  it  can  generally  be  misleading  since  it  assumes  that 
the  instantaneous  frequency  of  the  light  is  equivalent  to  a  time-averaged  frequency. 
Nevertheless,  the  model  is  instructive  as  it  suggests  that  functions  with  larger  slopes 
will  have  better  conversion  efficiency  from  FM  to  AM,  and  that  a  function  with  many 
large  high  order  derivatives  will  distort  the  AM  signal  more  than  one  with  a  more 
“linear”  function. 

The  optical  filter  used  as  an  FM  discriminator  impacts  the  performance  of  the 
FM-DD  link.  The  gain  of  the  fink  is  affected  by  the  discriminator’s  conversion  effi¬ 
ciency  from  FM  to  AM.  The  distortion  of  the  fink  is  limited  by  the  linearity  of  the 
conversion  process.  Because  of  the  long  coherence  lengths  of  lasers  used  in  communi- 
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cation  systems,  practical  filters  are  coherent  and  are  appropriately  analyzed  as  filters 
of  light’s  electric  field  amplitude  and  phase,  rather  than  as  filters  of  the  light’s  inten¬ 
sity.  In  this  work,  I  derive  guidelines  for  designing  discriminator  filters  to  ensure  low 
signal  distortion.  A  filter  with  perfectly  linear  phase  and  whose  electric  field  trans¬ 
mission  ramps  linearly  with  frequency  will  linearly  convert  FM  to  AM.  Because  the 
signal  is  AM  rather  than  IM,  the  photodetection  process  will  create  second-harmonics 
of  the  signal’s  Fourier-frequency  components.  This  can  be  dealt  with  by  a  balancing 
approach.  I  also  find  that  directly  converting  from  FM  to  IM  cannot  be  performed 
in  a  linear  fashion. 


1.4  Development  of  FM  discriminators 

Mach-Zehnder  interferometers  (MZI)  were  first  suggested  as  FM  discriminators,  but 
they  create  large  amounts  of  signal  distortion.  Alternatives  to  MZIs  have  been  pro¬ 
posed,  including  Fabry-Perot  interferometers,  fiber  Bragg  gratings  (FBGs)  and  ring 
resonators. 

The  work  of  Harris,  [11],  was  the  earliest  use  of  a  quadrature  biased  Mach  Zehnder 
interferometer  structure  to  discriminate  optical  FM.  An  interferometric  path  differ¬ 
ence  was  created  by  passing  the  light  through  a  birefringent  crystal  when  the  light’s 
polarization  was  angled  between  the  fast  and  slow  axes  of  the  crystal.  It  was  noted 
that  optimal  FM  to  AM  conversion  occurs  at  the  quadrature  bias  point.  The  technique 
was  also  applied  to  phase  modulated  light  in  [10].  Besides  PM  to  AM  discrimination, 
suppression  of  unwanted  incident  AM  was  done  by  applying  a  180  degree  phase  shift  to 
one  of  the  two  complementary  polarization  states  at  the  output  of  the  discriminator. 
The  initial  AM  canceled  when  both  polarization  states,  now  with  their  PM  in  phase 
but  AM  180  degrees  out  of  phase,  were  detected  at  a  single  polarization-insensitive 
photodetector. 
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The  authors  understood  that  linear  demodulation,  required  for  high  fidelity  signal 
transmission,  could  be  accomplished  with  a  discriminator  that  has  a  linear  FM  to 
AM  transfer  function,  and  that  high-order  filters  could  be  used  to  implement  this 
linear-field  discriminator.  In  [12],  they  proposed  a  linear-field  discriminator  using  a 
network  of  birefringent  crystals.  The  device  was  a  tenth-order  finite-impulse-response 
(FIR)  filter.  The  series  of  crystals  worked  as  a  series  of  cascaded  Mach  Zehnder  inter¬ 
ferometers,  and  the  network  was  equivalent  to  a  lattice  filter  architecture.  The  filter 
coefficients  chosen  were  the  exponential  Fourier  series  approximation  to  a  triangular 
wave. 

Another  physical  implementation  of  the  MZI  style  discriminator  using  mirrors  and 
beam  splitters  was  suggested  by  [13].  In  this  case,  balanced  photodetection  was  used 
to  cancel  AM.  Such  an  interferometer  was  experimentally  verified  by  [14].  [13]  also 
suggested  the  use  of  balanced  detection  for  the  birefringent  crystal  device  of  [11]. 

Concurrent  to  the  development  of  direct  frequency  modulation  of  semiconductor 
lasers  in  works  such  as  [15],  [16]  performed  digital  data  transmission  experiments  using 
a  Michelson  interferometer  to  discriminate  optical  frequency  shift  keying  (FSK). 

The  use  of  FM  semiconductor  lasers  and  discriminator  detection  was  extended  to 
transmitting  subcarrier-multiplexed,  analog  signals  for  applications  in  cable  television 
distribution.  Experimental  results  for  a  Fabry-Perot  discriminated,  FM  subcarrier- 
multiplexed  system  were  presented  by  [17].  An  array  of  optical  frequency  modulated 
DFB  lasers  and  a  Fabry-Perot  discriminator  were  used  to  transmit  and  demodulate 
a  large  number  of  microwave  FM,  analog  video  channels.  A  similar  system  was  also 
used  to  transmit  subcarrier- multiplexed,  digital  signals  in  [18]. 

Because  analog  links  require  high  linearity  and  low  noise,  a  number  of  authors  have 
derived  figures  of  merit  for  the  performance  of  analog  FM-DD  finks.  [19]  analyzed  the 
frequency-dependent  response  of  a  link  with  a  quadrature  biased  MZI  discriminator 
subject  to  large  modulation-depth  AM  and  FM.  [20]  studied  the  intermodulation 
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distortion  for  a  Fabry-Perot  discriminated  link  with  a  large  number  of  channels,  while 
taking  into  account  both  FM  and  IM  on  each  channel.  [21]  derived  figures  of  merit 
for  the  dynamic  range  of  a  phase  modulated  link  with  an  MZI  discriminator  and 
balanced  detection. 

[22]  studied  a  fink  with  an  arbitrary  discriminator.  The  general  formulae  were 
applied  to  the  particular  cases  of  an  MZI  and  a  Fabry-Perot  interferometer.  However, 
the  analysis  was  inaccurate  since  it  looked  at  the  system  in  terms  of  fight  intensity 
transmission  through  the  interferometer.  The  transmission  was  expanded  in  terms 
of  a  Taylor  series.  The  analysis  assumed  that  derivatives  of  the  transmission  spec¬ 
trum  of  the  interferometer  (in  the  Fourier-frequency  domain)  with  respect  to  the 
instantaneous  optical  frequency  were  proportional  to  overall  link  nonlinearity.  Sim¬ 
ilar  (inaccurate)  theoretical  analyses  using  Taylor  series  were  published  by  [23]  and 
[24].  However,  these  papers  did  include  new  models  for  the  nonlinearities  in  the  lasers’ 
FM  and  included  the  effects  of  residual  IM. 

To  improve  the  linearity  of  an  FM-DD  fink,  many  alternatives  to  the  Mach- 
Zehnder  and  Fabry-Perot  interferometers  have  been  suggested.  Except  for  the  early 
work  of  [12],  all  of  these  “linearized”  discriminators  were  designed  such  that  the  fil¬ 
ter’s  optical  intensity  transmission  was  linear  with  frequency,  rather  than  the  field 
amplitude.  [25]  and  [26]  proposed  pairs  of  chirped  fiber-Bragg  gratings  with  either 
the  index  variation  or  chirp  rate  varied  nonlinearly.  [27]  proposed  a  frequency  dis¬ 
criminator  based  on  an  MZI  with  ring  resonators  in  its  arms.  [28]  suggested  that  the 
linearity  of  a  Sagnac  discriminator  could  be  improved  by  adding  ring  resonators. 

Experimental  and  theoretical  results  using  fiber-Bragg  gratings  were  presented  in 
[29,  30,  31,  32,  33,  34].  These  experiments  used  pairs  of  complementary  gratings 
designed  to  have  linear  intensity.  The  gratings  were  low-biased  to  perform  carrier 
suppression.  In  [31,  33],  the  authors  presented  a  clipping-free  dynamic  range  limit 
for  an  EM-DD  system.  (In  related  work,  [35,  36],  the  authors  used  Bragg  gratings 
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to  convert  phase  modulation  into  single  sideband  modulation.)  The  authors  later 
realized  that  linear  FM  to  IM  discrimination  was  not  consistent  with  their  theoretical 
analysis,  [32]: 

[...]  the  ideal  linear  power  reflectivity-versus-frequency  curve  does  not 
result  in  an  ideal  half-wave  rectification,  as  suggested  by  the  simple 
time-domain  view.  Rather,  in  addition  to  the  signal  component,  the 
output  includes  a  dc  component  as  well  as  a  nonlinear  distortion. 

They  explained  the  discrepancy,  [34]: 

The  reason  this  intuition  fails  is  that  combining  a  time-domain  view  of  the 
FM  signal  (instantaneous  frequency,  not  averaged  over  time)  with  a 
frequency  domain  view  of  the  FBG  filter  response  is  inconsistent  with 
the  frequency  domain  analysis  [...] 

It  is  erroneous  to  think  of  the  modulated  signal  in  terms  of  its  instantaneous  frequency 
while  looking  at  the  frequency  spectrum  of  the  filter.  The  carrier  is  not  really  being 
swept  along  the  ramp  of  the  filter  by  the  modulation,  so  analyzing  it  in  the  same  way 
as,  for  example,  the  small-signal  current  to  voltage  relationship  of  an  amplifier  is  not 
correct.  In  this  work,  I  present  FM  to  AM  discrimination  as  an  alternative  which 
can  produce  a  microwave  photonic  link  with  low  distortion.  The  following  chapter 
analytically  demonstrates  this  result. 
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Chapter  2 


Analytical  link  analysis 


In  this  chapter,  I  derive  figures-of-merit  for  an  FM-DD  link  that  uses  an  arbitrary  op¬ 
tical  filter  for  discrimination  [37].  I  find  expressions  for  the  currents  at  each  microwave 
frequency  at  the  output  of  the  link  under  a  two-tone  test.  I  take  a  small-modulation- 
depth  approximation  and  obtain  expressions  for  the  signal-to-noise  ration  (SNR), 
second-order  and  third-order  output  intercept  points  (OIP2  and  OIP3),  spurious-free 
dynamic  range  (SFDR)  and  noise  figure  (NF).  I  apply  these  general  formulae  to  the 
specific  cases  of  the  Mach  Zehnder  interferometer,  a  linear-intensity  filter  and  a  linear- 
field  filter.  The  analysis  of  the  MZI  is  consistent  with  earlier  theory.  I  show  that  the 
linear-intensity  filter  does  not  convert  FM  to  IM  in  a  linear  fashion,  but  the  linear- 
field  filter  actually  does  convert  FM  to  AM  in  a  linear  fashion.  For  the  linear-field 
filter,  I  derive  the  noise  figure’s  dependence  on  the  link’s  regime  of  operation  and 
quantify  the  effect  of  the  laser’s  residual  IM  on  the  distortion. 

2.1  Response  of  filtered  FM  link 

In  this  section,  I  derive  expressions  for  the  rf  photocurrents  produced  at  the  out¬ 
put  of  an  FM-DD  fink  with  an  arbitrary  filter.  An  FM-DD  fink  consists  of  an  FM 
laser,  an  optical  filter  and  a  photodetector.  The  source’s  residual  IM  and  the  filter’s 
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nonlinearities  are  the  main  constributions  to  distortion  on  the  link.  To  quantify  the 
nonlinearities,  I  apply  a  two-tone  distortion  test.  An  optical  signal  that  is  phase  or 
frequency  modulated  by  two  sinusoidal  tones  can  be  represented  by  the  time  varying 
electric  field 


e  (t)  =  cos  [2TTfct  +  (3i  sin  (27r/it)  -f  sin  (27r/2t)]  (2.1.1) 

where  Popt  is  the  rms  optical  power,  k  is  a  constant  relating  optical  field  and  optical 
power  such  that  Popt  =  (e  (t)^)  / fc  is  the  frequency  of  the  optical  carrier,  /i  and 
/2  are  the  modulation  frequencies  and  (3i  and  (32  are  the  angle  modulation  depths. 
For  FM,  each  modulation  depth  is  equal  to  the  maximum  optical  frequency  deviation 
of  the  carrier  induced  by  the  modulation  divided  by  the  frequency  of  the  modulation, 
P  =  5f/  f  .  The  modulation  of  the  light  can  be  thought  of  in  terms  of  variations  in  the 
instantaneous  frequency  of  the  light  due  to  the  applied  signal.  The  optical  frequency, 
or  wavelength,  varies  sinusoidally  in  time.  The  instantaneous  frequency  of  the  light 
is  given  by  the  derivative  of  the  phase  of  the  fight. 


[271  fct  +  Pi  sin  {271  fit)  +  P2  sin  (27r/2t)]  =  /c  +  5/i  cos  (27r/it)  +  5 cos  {27if2t) 

l7i  at 

(2,1,2) 

An  FM  laser  is  non-ideal  as  it  produces  undesired  residual  IM  and  includes  noise. 
The  correction  to  the  electric  field  is 

e  {t)  =a  {t)  +  KiJ 2Popt  [l  +  n  (t)]  (2.1.3) 

■  \/l  +  mi  cos  {271  fit  +  P)  +  1712  cos  {277 f2t  +  p) 

■  cos  [27r/ct  -|-  Pi  sin  (27r/it)  -f  P2  sin  {27if2t)  +  ^{t)] 

where  n  {t)  is  the  RIN  of  the  source,  (p{t)  is  the  phase  noise  of  the  source,  a  {t)  is  the 
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ASE  noise  from  an  optical  amplifier,  mi  and  m2  represent  tfie  IM  deptfis  for  tfie  two 
tones  and  0  is  tfie  pfiase  difference  between  tfie  IM  and  tfie  FM.  Tfie  fink  will  also 
amplify  tfiermal  noise  present  at  tfie  input. 

An  arbitrary  optical  filter  is  used  on  tfie  fink  to  convert  FM  to  IM.  Using  tfie 
Jacobi-Anger  expansion,  the  electric  field  after  the  filter  can  be  expressed  as  an  in- 
finited  weighted  sum  over  sidebands.  I  employ  a  shorthand  notation  to  describe  the 
electric  field  transmission  at  each  frequency  in  the  optical  spectrum  that  corresponds 
to  an  optical  sideband, 

hn,p  =  h{f  =  fc  +  nfi+  p/2)  (2.1.4) 

where  n  and  p  are  integer  indices  and  h  is  the  complex  transfer  function  of  the  filter, 
representing  its  phase  and  amplitude  response,  including  any  insertion  losses  or  optical 
amplifier  gain.  For  example,  /ro,o  is  the  field  transmission  for  the  optical  carrier,  and 
/r_i  o  is  the  transmission  of  the  negative,  first  order  sideband  spaced  fi  away  from  the 
carrier. 

After  the  filter,  the  fight  is  incident  upon  a  photodetector.  I  derive  the  output 
current  from  the  photodetector  and  approximate  it  for  small  modulation  depth.  The 
standard  definitions  for  the  linearity  figures  of  merit  rely  on  this  approximation.  The 
residual  IM  depth  and  the  intensity  noise  are  assumed  to  be  much  smaller  than  the 
FM,  so  the  square  root  in  (2.1.3)  can  be  expanded  using  a  Taylor  series,  yielding 


e  (t)  (t)  +  K^/2Popt 


(2.1.5) 


111 
1  +  -mi  cos  (271  fit  +  0)  +  -m2  cos  (271  f2t  +  0)  +  -n  (t) 


cos  [2n fct  +  01  sin  (27r/it)  +  02  sin  (27rf2t)  +  (p(t)] 


Ignoring  noise,  this  can  be  written  using  an  angular  addition  trigonometric  identity 
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as 


e  (t)  =  K^J2PoptRe\^ 

cos  [2nfct  +  (3i  sin  (27r/it)  +  (32  sin  (27r/2t)] 

+  cos  [27r  (/c  +  /i)  t  +  A  sin  {2nfit)  +  sin  {27rf2t)  +  0] 

+  cos  [27r  (fc  -  fi)t  +  /3i  sin  (27r/it)  +  /^a  sin  (27r/2t)  -  0] 

+  ^1712  cos  [27r  (/c  +  /2)  t  +  I3i  sin  {2nfit)  +  /^a  sin  (27r/2t)  +  0] 

+  ^7712  cos  [27r  (/c  -  /i)  t  +  /3i  sin  {2nfit)  +  /^a  sin  (277/2^)  -  0]  | 

The  Jacobi-Anger  expansion  is  given  by  =  Yl'^=-oo3^'^n  (z)  where  j  is  the 

imaginary  unit  and  Jn{z)  is  a  Bessel  function  of  the  first  kind.  Applying  this  formula, 
the  expression  expands  to 


e  (t)  =  K^j2PoptRe< 

00  00 

E  E  Jn  (A)  Jp  W2)  exp  [j27r  {fc  +  nfi  +  p/a)  t] 

:  — O 
1 


n=—oo  p=—oo 

00  00 


^  ^  (^1)  '^p  (^2)  exp  [j27r  {fc  +  [n  +  l]/i  +  p/2)  t  +  //] 


n=— 00  p=—oo 
00  00 


+  -mi  E  E  Jn  (A)  Jp  (A)  exp  [j27r  {fc  +  [n  -  l]/i  +  p/a)  t  -  //] 


n=—oo  p=—oo 
00  00 


^272^^  (^2)  exp  [j27r  {fc  +  nfi  +  [p  +  Ij/a)  t  +  //] 


n=— 00  p=—OQ 
00  00 


+  7"^2  X]  X]  A(A)  A(A)exp[j27r(/c  +  n/i  +  [p- l]/2)t- j0] 


n=— 00  p=—oo 


The  signal  passes  through  the  FM  discriminator  filter.  The  electric  field  after  the 
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filter  is 


e  (t)  =  K^j2PoptRe< 

OO  OO 

E  E  Jn  (A)  Jp  {(^2)  hn,p  exp  [j2TT  (/c  +  ufi  +  p/2)  t] 


n=—oo  p=—oo 

CXD  00 


+  -mi  E  E  Jn  (A)  Jp  {1^2)  fin+l,p exp  [j2Tl  (/c  +  [n  +  l]/i  +  p/2)  t  +  j0] 


n=—oo  p=—oo 

OO  OO 


+  -mi  E  E  Jn  (/3l)  (/p  ^1^2)  ^n— l,p^^P  [j27r  (/c  “h  [xi  l]/l  “1“  P/2)  t  j4^] 


n=—oo  p=—oo 

OO  OO 


+  -m2  E  E  Jn  (/3l)  Jp  {/32)  hn,p+l  exp  [j27l  (/c  +  n/i  +  [p  +  l]/2)  t  +  //] 


n=—oo  p=—oo 

OO  OO 


+  -m2  E  E  Jn  (A)  (/52)  fin,p-l  exp  [j27r  (/c  +  ufi  +  [p  -  l]/2)  t  -  j0] 


n=— OO  p=—oo 


The  indices  of  each  infinite  sum  can  be  renumbered  to  obtain 


e  (t)  =  K^J2PoptRe< 

OO  OO 

E  E  Jn  (A)  Jp  {(^2)  hn,p  exp  [/27r  (/c  +  nfi  +  p/2)  t] 


n=—oo  p=—oo 

-  0000 


+  -mi  E  E  Jn-i  iPi)  Jp  {P2)  hn,p  exp  [j27r  (/c  +  n/i  +  p/2)  t  +  j/] 


n=—oo  p=—oo 

OO  OO 


+  -m,  E  E  Jn+1  iPi)  Jp  {132)  hn,p  exp  [/27r  {fc  +  nfi  +  p/2)  t  -  //] 


n=—oo  p=—oo 

OO  OO 


+  -m2  E  E  Jn  {(3i)  Jp-I  {(32)  hn,p  exp  [/27r  {fc  +  nfi  +  p/2)  t  +  jf] 


n=—oo  p=—oo 

OO  OO 


+  -m2  E  E  Jn  {Pi)  Jp+I  {P2)  hn,p  exp  [/27r  {fc  +  nfi  +  p/2)  t  -  jf] 


n=—oo  p=—oo 


This  simplifies  to  a  compact  expression  for  the  signal  after  the  filter  in  terms  of  its 
frequency  components, 
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5(t)  =  Ky^W^tRe 


oo 


oo 


2^  2^ 

,  n=—oo  p=—oo 


jn,p  exp  [j27r  (/c  +  nfi  +  p/2)  t]  i  (2.1.6) 


where  I  define 


jn,p=hn,p{Jn{f3l)Jp{/32)  (2-1-7) 

+  [Jn-1  iPl)  +  Jn+1  (/3l)  e  Jp  {(32) 

+  -^^2Jn  (/^i)  ['7p-i  (/52)  +  '7p+i  (/3i)  e  j- 


The  electric  field  is  incident  upon  a  photodetector  at  the  termination  of  a  fiber-optic 
fink.  The  photodetector  is  assumed  to  be  an  ideal  square-law  detector  operating  in 
its  linear  regime  with  responsivity  3?  .  The  photocurrent  is 


i{t)  =  ^Popt  j  jn,pjlk  exp  [j27r  {[n  -  g]  fi  +  \p  -  k]  /a)  t]  I 

\^n=—oop=—oog=—ook=—oo  J 

(2.1.8) 

This  can  be  split  up  into  the  dc  term,  harmonics  of  /i,  harmonics  of  /2  and  mixtures 
between  fi  and  /2. 


OO  OO 


i{t)  =^Popt  < 

\g=—(yD  k=—oo 
oo,n^g  00  00 

+  E  E  Z  g,k  exp  [j27r  [n  -  g]  fp] 

n=—(yo  g=—oo  k=—oo 

oo,p^k  00  00 

+  E  E  E  jg,pj*g,k  exp  [j27r  \p  -  k]  f2t] 

p=—oo  g=—oo  k=—oo 


oo,n^g  oo,p^k  oo  oo  ^ 

E  E  E  E  Jn,pjlk  exp  [j27r  {[n-g]fi  +  [p-  k]  /2)  t]  I 

n=—oop=—oog=—ook=—oo  J 
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The  indices  of  each  infinite  sum  can  be  renumbered  to  obtain 


oo  oo 


Y.  Y  1^9 ‘I" 

\^g=—oo  k=—oo 
oo,n^0  oo  oo 

+  E  E  E  Jn+g,kj*g,k  exp  [j27rn/it] 

n=— OO  g=—OQ  k=—oo 
oo,p7^0  oo  oo 

+  E  E  E  jg,P+kjg,k  exp  [j2'Kpf2t] 


p=—OQ  g=—OQ  k=—oo 
oo,n7^0  oo,p7^0  oo 


+  E  E  E 


jn+g,p+kjlk  exp  [j27r  (n/i  +  p/2)  t] 


n=—oo  p=—oo  g=—oo  k=—oo  ) 

The  double  infinite  sums  over  n  and  p  are  rewritten  as  singly  infinite  sums,  and  the 
sums  over  negative  integers  have  their  signs  Hipped  giving 


\^g=—oo  k=—oo 
00  00  00 

+EE  E  {jn+g,kjlk  exp  [j27rn/it]  +  j-n+g,kj*g,k  exp  [-j27rn/it]) 

n=l  g=—oo  k=—oo 
00  00  00 

+EE  E  {jg,p+kj*g,k  exp  [j27rp/2t]  +  jg,-p+kjlk  exp  [-j27rp/2t]) 

p=l  g=—oo  k=—oo 
00  00  00  00 

-^EEE  E 

n=l  p=l  g=—oo  k=—oo 

{jn+g,p+kjlk  exp  [j27r  (n/i  +  p/2)  t]  +  j-n+g,-p+kj*g,k  exp  [-/27r  (n/i  +  p/2)  t] 
+jn+g,-p+kj*g^k  exp  [/27r  (n/i  -  p/2)  t]  +  j_n+g,p+kjlk  exp  [-/27r  (n/i  -  p/2)  t])  } 

A  number  added  to  its  complex  conjugate  is  twice  the  real  part.  With  this  simpli¬ 
fication,  this  arranges  to  a  final  expression  for  the  fink  output  given  an  arbitrary 
filter: 
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i(t)  =‘RPoptRe 


OO 


2 


(2.1.9) 


OO 

\h,k 

g=—OQ  k=—oo 

OO  OO  OO 

+2E5:  E  jn+g,kjlk  exp  [j27rn/it] 

n=l  g=—oo  k=—oo 

OO  OO  OO 

+  2EE  E  jg,p+kj*g,k  exp  [j2TTpf2t] 

p=l  g=—oo  k=—oo 

OO  OO  OO  OO 

+  2EEE  E  jn+g,p+kfg,k  exp  [j2TT  {ufi  +  p/a)  t] 

n=l  p=l  g=—oo  k=—oo 

OO  OO  OO  OO 

+2  EEE  E  jn+9,-p+kj*g,k  exp  [i27r  (n/i  -  p/a)  t] 

n=l  p=l  g=—oQ  k=—oo 

The  double-sum  over  indices  g  and  k  gives  the  contribution  of  each  pair  of  optical 
sidebands  that  beat  together  to  produce  the  rf  photocurrent.  In  this  form,  the  current 
is  separated  into  different  frequency  components  which  are  indicated  by  the  summa¬ 
tion  indices  n  and  p  .  The  first  term,  where  n  and  p  are  both  identically  zero,  gives  the 
dc.  The  second  term,  a  summation  over  the  index  n,  gives  the  fundamental  tone  at 
frequency  f\  and  its  harmonics.  The  third  term,  a  summation  over  the  index  p,  gives 
the  fundamental  tone  at  frequency  /2  and  its  harmonics.  The  fourth  term  is  the  sum 
frequencies  produced  by  the  mixing,  and  the  fifth  term  is  the  difference  frequencies 
produced  by  the  mixing. 

For  small  modulation  depth,  /3  1,  and  no  residual  IM,  m  =  0,  the  Bessel 

functions  can  be  approximated  by  Jo{z)  ~  1  and  Jn{z)  (;2/2)'”'  /  |n|!,  for  positive 
n,  noting  that  J-n{z)  =  (—1)”  Jn{.z).  Keeping  terms  of  lowest  polynomial  order,  the 
current  simplifies  to  the  following  equation  (2.1.10).  This  equation  gives  the  small 
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signal  approximation  for  any  frequency: 


i{t)  =3f?Popii?e  {|/ro,o|^ 

/sr  (-1)" 


(2.1.10) 


^=i  g=i 
oo  J)  r,p 


2EE 


(-i)‘ 

2P  {p  —  k)\k\ 


ho^p_khl_k  exp  [j2TTpf2t] 


p=l  fc=0 

oo  oo  n  p  on  oP  / _ -i  \9+^ 

2  E  E  E  E  “p  12^’^  ‘1 

n=l  p=l  p=0  k=0  ^  n  J 


^=i  p=i  g=ij  k=' 
oo  oo  n  P  on  op 


+2EEEE 

n=l  p=l  p=0  k=0 


(-1) 


p+p+fc 


2^+P  {n-g)\g\  {p  -  k)  k\ 


hn-g-p+kh*_  k  exp  [j27r  (n/i  -  p/2)  t] 


There  are  four  current  components  of  interest.  The  amplitude  of  the  dc,  as  should 
be  expected,  is  proportional  to  the  optical  power  in  the  optical  carrier  after  the 
filter.  The  current  at  the  fundamental  frequency  /i  is  linearly  proportional  to  the 
modulation  depth.  It  depends  on  the  negative  and  positive  first-order  sidebands 
beating  with  the  optical  carrier.  The  current  at  the  second-harmonic  frequency  2/i 
has  a  quadratic  relationship  to  modulation  depth.  It  depends  on  the  second-order 
sidebands  beating  with  the  optical  carrier,  as  well  as  the  first-order  sidebands  beating 
with  each  other.  The  current  produced  at  the  difference  frequency  2/i  —  /2  is  a 
third-order  intermodulation  product.  These  currents  are 


idc=^PoptXo  (2.1.11) 

ip  =^Poptl3iRe  {Xi  exp  [/27r/it]}  (2.1.12) 

*2/1  =^Popt^l3iRe  {X2  exp  [j47r/it]}  (2.1.13) 

*2/1-/,  =^Poptlf3^^f32Re{Xsexp[j27r{2h-  f2)t]}  (2.1.14) 
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where  for  convenience,  I  define  the  following  complex  constants 


-^0  =^0, 0^0,0 

(2.1.15) 

0^0,0  ~ 

(2.1.16) 

=fii, 0^0,0  + 

(2.1.17) 

X2  =^2, 0^0,0  ~  +  ho,oh-2,0 

(2.1.18) 

X3  =  —  h2-ih*Q  Q  +  ^2, 0^0,1  A  Q 

(2.1.19) 

+  fio,ofil2,i  ~  hQ-ih*_2  Q  —  i 

Each  rf  photocurrent  outputs  an  rms  power  into  the  load  resistance,  Rioad,  that  is 
proportional  to  the  square  of  the  dc  current. 


P/,  lA'il" 

(2.1.20) 

P2/.  |A'2|= 

(2.1.21) 

2/1 -/2 

(2.1.22) 

In  this  section,  I  have  derived  closed  form  expressions  for  the  photocurrents  at 
different  frequencies  at  the  output  of  a  filtered  FM  link.  A  general  result  has  been 
given  in  (2.1.9)  which  includes  residual  intensity  modulation,  and  can  be  solved  to 
arbitrary  precision  by  taking  a  large  number  of  terms  in  the  infinite  sum.  A  small 
signal  approximation,  (2.1.10),  gives  the  output  current  at  any  frequency  component 
of  interest.  Expressions  for  the  photocurrent  at  the  fundamental,  second  harmonic 
and  third  order  intermodulation  distortion  have  been  derived,  which  will  be  useful  in 
expressing  figures  of  merit  for  distortion  and  dynamic  range. 
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2.2  Signal  to  noise  ratio 

In  this  section,  I  derive  the  signal  to  noise  ratio  (SNR)  for  the  small  signal  approxi¬ 
mation  of  an  arbitrary  link.  A  passive  link  with  no  amplification  will  be  considered, 
so  the  primary  noises  seen  at  the  detector  are  shot,  thermal,  phase  and  RIN.  The 
shot  noise  spectral  density  is  proportional  to  the  dc  from  the  photodetector  and  q, 
the  elementary  charge.  The  thermal  noise  spectral  density  is  equal  to  the  product 
Boltzmann’s  constant,  ks,  and  the  temperature,  Tk- 

Sgn  ‘^Q^dcRload 

Stn  =kBTK  (2.2.2) 


Assuming  a  Lorentzian  model  for  the  laser’s  spectral  line,  the  phase  noise  on  the 
optical  carrier  is  white  noise  with  spectral  density  proportional  to  the  laser’s  3-dB 
linewidth,  Au  [38].  The  phase  fiuctuations  are  converted  to  intensity  fiuctuations 
by  the  filter  in  the  same  manner  as  it  converts  the  modulation.  The  average  phase 
fiuctuations  in  a  small  bandwidth  near  some  frequency,  /,  are 


{‘fitf) 
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Near  the  first  modulation  frequency,  /i,  the  power  spectral  density  of  the  phase  noise 
is 

Spn  ~  Rload^^P^pt^  |Xi|'  =  (2.2.3) 

TT/i  TT/i 

The  modulation  is  assumed  to  be  below  the  relaxation  frequency  of  the  laser,  so  the 
RIN  is  modeled  as  white  noise.  The  power  spectral  density  of  the  noise  at  the  output, 
near  the  modulation  frequency  is 


~  ^RloadAc^O 


-2 


B 


(2.2.4) 
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where  B  is  the  bandwidth  in  Hz.  The  total  noise  power  is 


^  noise  ^ 


{Ssn  +  Stn  +  Spn  +  Sin)  B 


(2.2.5) 


The  signal  to  noise  ratio  (SNR)  is  Pf^/Pnoise-  If  the  SNR  is  phase  noise-limited, 
which  is  the  case  for  sufficient  optical  powers,  moderate  RIN,  and  efficient  conversion 
of  phase  fluctuations  into  intensity  fluctuations,  the  SNR  is  given  by 

SNR  =  h^fTi/AuB  (2.2.6) 

For  FM  discriminators  that  have  constant  conversion  efficiency  over  a  large  enough 
bandwidth,  the  upper  bound  on  the  SNR  for  an  FM-DD  link  is  determined  by  the 
modulation  depth  and  linewidth  of  the  FM  laser. 


2.3  Distortion 

The  signal  distortion  caused  by  the  FM-DD  link  can  be  described  by  the  output  power 
at  frequencies  that  are  harmonics  and  mixing  terms  of  the  modulation  frequencies. 
For  now,  I  assume  there  is  no  residual  IM.  If  each  of  the  two  modulation  tones  has 
equal  modulation  depth,  the  IMD’s  power,  (2.1.22),  is  equal  to  the 

signal  power,  (2.1.20),  for  modulation  depth  5j  =  8/1/2  |-^i|  /  The  corresponding 
third-order  output  intercept  point  (OIP3)  is 

OIP3  =  4Rw*LXo-2/rV2  |Xi|=*  /  IX3I  (2.3.1) 

The  second  harmonic’s  power,  (2.1.21),  is  equal  to  the  signal  power,  (2.1.20),  for 
modulation  depth  Sj  =  16fi  \Xif  /  \X2f .  The  corresponding  second-order  output 
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intercept  point  (OIP2)  is 


OIP2  =  SRioad^dcXo^  |Xi|"  /  IX2I'  (2.3.2) 

For  an  arbitrary  filter,  the  distortion  will  depend  on  the  particular  modulation  fre¬ 
quencies  chosen.  One  desires  to  maximize  |Xi|  and  minimize  IX2I  and  IX3I  to  reduce 
the  distortion.  A  link  with  zero  IX2I  or  IX3I  will  have  infinite  OIP2  or  OIP3. 


2.4  Spurious  free  dynamic  range 


The  spurious  free  dynamic  range  (SFDR)  is  defined  as  the  SNR  at  the  maximum 
usable  modulation  depth.  This  can  be  defined  when  either  the  second-order  or  third- 
order  distortion  products  breach  the  noise  fioor.  For  a  phase  noise-limited  link,  the 
IM3  is  equal  to  the  noise  power  at  modulation  depth  /  IX3 

For  a  link  limited  by  IM3,  using  (2.2.6),  the  SFDR  is 


SFDR^ 


t  B/\v  lA'al  ) 


(2.4.1) 


For  a  phase  noise-limited  link,  the  power  at  the  second-harmonic  frequency  is  equal 
to  the  noise  power  at  modulation  depth  =  4/i |Xi|  /  IX2I.  For  a  link  limited 
by  the  second-harmonic  distortion,  the  SFDR  is 


SFDR2 


2/iv/^|Xi 

v/ATRiXsl 


(2.4.2) 


These  figures-of-merit  are  often  defined  with  respect  to  1  Hz  bandwidth.  They  gener¬ 
ally  depend  on  the  particular  modulation  frequencies  chosen.  Maximizing  the  ratios 
of  |Xi|  /  1X2 1  and  |Xi|  /  IX3I  will  improve  the  dynamic  range  of  the  link. 


2/3 
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2.5  Specific  FM  discriminators 


Because  the  bandwidth  of  optical  systems  is  very  large,  the  spectrum  of  optical  filters 
can  include  many  repeating  periods.  To  specify  the  filter  response,  one  starts  with  the 
frequency  response  of  the  filter  normalized  over  one  free  spectral  range  (FSR),  written 
in  terms  of  angular  frequency  =  — tt  to  tt.  One  period  of  the  filter  is  centered  at  a 
chosen  center  frequency  /o,  which  is  not  necessarily  the  same  as  the  carrier  frequency 
fc-  I  define  the  bias-frequency  offset  of  the  filter  as  fb  =  fc  —  fo,  which  can  be  adjusted 
by  tuning  the  wavelength  of  the  transmitting  laser  or  the  physical  parameters  of  the 
optical  filter.  The  amplitude  of  the  transfer  function  corresponds  to  the  electric  field 
transmission  of  the  filter.  The  maximum  bounds  on  the  amplitude  are  -1  and  1. 
Negative  transmission  corresponds  to  a  phase  shift  at  zero  field  transmission.  The 
phase  shifts  between  adjacent  periods  of  the  filter  are  determined  by  the  filter  type. 

2.5.1  Mach-Zehnder  interferometer 

The  simplest  filter  used  as  an  FM  discriminator  is  an  MZI  with  50%  coupling  ratios. 
One  arm  of  the  interferometer  has  a  time  shift  with  respect  to  the  second  arm.  A 
normalized  period  of  the  filter  is  written  as 

h{^)  =  ^-^exp{-jQ)  (2.5.1) 

The  filter  is  typically  biased  at  quadrature,  giving  in  our  notation  (see  equ.  2.1.4)  a 
transfer  function  of 


hn,p  =  ^  ^  (n/i  -f  p/2)  r]  (2.5.2) 

where  r  is  the  time  delay  between  the  two  arms.  (This  is  not  written  as  h{f)  in 
order  to  avoid  confusion  between  modulation  and  optical  frequencies).  The  intensity 
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response  is  a  sinusoid 


hn,pK^p  =  -{1- sin  [-2%  (n/i  +  p/2)  r]} 


(2.5.3) 


Using  the  transfer  function,  I  evaluate  the  link  constants 


Xn  =:r 


(2.5.4) 


X2  =0 

X3  =  -  (1  -  (1  - 


(2.5.5) 

(2.5.6) 

(2.5.7) 


As  expected  for  an  MZI  at  quadrature,  I  find  that  there  is  no  second-harmonic  so 
that  OIP2  is  infinite.  Using  the  approximation  /ir, /2r  1,  the  absolute  value  of 

the  other  coefficients  are 


|Xi|  =sin  (tt/it)  tt/iT 


1X3 1  =4sin^  (tt/it)  sin  (7r/2r)  4:TT^fif2r^ 


(2.5.8) 

(2.5.9) 


The  power  at  the  fundamental  and  IMD3  frequencies  are 


P/,  =2Rioad'i^dc^%'^'^r^ 


P2h-h  =7,Ploadi^dJ%S%TT^T^ 


(2.5.10) 

(2.5.11) 


The  OIP3  of  the  quadrature  biased  MZI  and  the  phase  noise-limited  SFDR  are 


OIP3  —4:Rioadidc 


SFDRs  =  (Az/Pttt^) 


2\-2/3 


(2.5.12) 

(2.5.13) 


27 


Both  the  useful  bandwidth  and  the  SFDR  are  improved  by  having  a  short  time  delay 
in  the  MZI.  Identical  results  to  (2.5.10)  and  (2.5.12)  are  found  by  [21]  (with  prefactors 
of  x4  due  to  twice  the  current  with  balanced  detection),  which  supports  the  general 
analysis. 

2.5.2  Linear  intensity 

A  number  of  groups  have  proposed  or  built  optical  filters  that  have  a  transfer  function 
linear  in  optical  intensity  versus  frequency  and  small  group  delay.  Within  one-half 
period,  the  transfer  function  can  be  represented  by 

hn,p  =  y^A{fb  +  nfi  +  p/2)  exp  [-j27r  (/b  -f  n/i  +  p/2)  r]  (2.5.14) 

where  A  is  a  slope  in  units  of  inverse  frequency  and  r  is  a  time  delay.  The  intensity 
response  is 

hn,pK^p  =  A{fb  +  nfi  +  p/2)  (2.5.15) 
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which  is  linear  in  slope  A.  Using  the  transfer  function,  I  evaluate  the  link  constants: 


(2.5.16) 

(2.5.17) 

(2.5.18) 

(2.5.19) 

(2.5.20) 


Generally,  X2  and  X3  are  non-zero  for  this  discriminator,  even  if  the  square  roots  are 
expanded.  This  means  that  an  FM  discriminator  that  is  linear  in  optical  intensity 
will  still  produce  second-order  and  third-order  distortion.  Mixing  in  the  photodetector 
produces  cross  terms  that  are  not  eliminated.  An  FM  discriminator  that  is  linear  is 
optical  intensity  will  not  produce  a  distortion-less  link. 


2.5.3  Linear  electric  field 

The  ideal  filter  for  an  FM-DD  link  is  an  optical  filter  that  is  linear  in  electric  field. 
I  find  that  this  filter  has  second  order-distortion  that  is  produced  by  the  square- 
law  detection,  but  no  other  higher-order  distortion.  Within  one  period,  the  field 
transmission  ramps  linearly  with  frequency,  and  the  filter  has  linear  phase.  The 
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transfer  function  is 


hn,p  =  A{fb  +  nfi  +  p/2)  exp  [-j27r  {fb  +  nfi  +  p/2)  r]  (2.5.21) 

where  A  is  a  slope  in  units  of  inverse  frequency  and  r  is  a  time  delay.  In  the  intensity 
domain,  the  filter  looks  quadratic. 


hn,pK^p  =  {fb  +  nfi  +  p/2)^  (2.5.22) 

The  fink  constants  are 

Xo  =A^f^  =  T  (2.5.23) 

Xi  =2A^  fifbC-^'^^^^^  =  2yl/iT^/2g-i2^/ir  (2.5.24) 

X2  (2.5.25) 

X3  =0  (2.5.26) 


I  define  the  constant  T  to  describe  the  dc  bias  of  the  filter,  which  is  the  fraction  of 
optical  power  transmitted  by  the  filter  at  the  optical  carrier  frequency.  I  find  the 
photocurrents  at  the  output  of  the  photodetector. 


idc  =^PoptT  (2.5.27) 

if^  =2idcT~^^^Sf^Acos  [2nfi  {t  -  r)]  (2.5.28) 

*2/i  =^idcT~^S\A'^  cos  [47r/i  {t  -  r)]  (2.5.29) 


The  magnitude  of  the  fundamental  current  is  linearly  proportional  to  the  slope  of 
the  filter  and  linearly  proportional  to  the  frequency  modulation  depth  6f  .  This 
filtered  link  produces  distortion  at  the  second-harmonic  of  each  modulation  tone. 
The  distortion  is  caused  by  the  first-order  sidebands  beating  with  each  other.  The  rf 
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power  at  the  fundamental  and  second-harmonic  frequencies  are 


Pt,  (2.6.30) 

P2/,  (2.6.31) 

In  the  small  modulation  depth  approximation,  this  ideal  FM-DD  link  has  no 
other  higher-order  distortion.  Using  a  symbolic  algebra  solver,  I  verified  that  the 
current  is  zero  for  all  intermodulation  and  harmonic  frequencies  up  to  sixth  order. 
At  a  given  harmonic,  sum  or  difference  frequency,  if  all  the  sidebands  in  the  sum  in 
(2.1.9)  corresponding  to  that  frequency  fall  within  a  region  of  the  filter  that  closely 
approximates  the  desired  linear  ramp  function,  the  output  current  is  zero. 

The  second-harmonic  can  be  suppressed  if  the  output  of  the  optical  system  is 
detected  using  balanced  detection.  The  filter  before  the  second  photodetector  is 
designed  to  have  a  slope  complementary  to  the  first  filter.  Its  transfer  function  is 

hn,p  =  A{fb-  nfi  -  p/2)  exp  [-j27r  {fb  +  nfi  +  p/2)  r]  (2.5.32) 

The  current  component  at  the  fundamental  frequency  will  be  180°  out  of  phase  be¬ 
tween  the  two  photodetectors,  but  the  second-harmonic  will  be  in  phase.  Subtracting 
the  second  current  from  the  first,  the  second-harmonic  will  cancel. 

Additional  sources  of  nonlinearity  are  the  frequency  modulated  laser  source,  op¬ 
tical  fibers  and  photodetector.  For  sufficient  modulation  depth,  the  dominant  FM 
sidebands  will  fall  outside  the  bandwidth  of  the  filter  and  this  saturation  will  cause 
nonlinearities. 
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2.6  Optimization  of  the  linear  link 

2.6.1  Gain,  noise  fignre  and  low  biasing 

Low  biasing  the  filter,  meaning  that  fc  is  very  close  to  the  zero  transmission  point 
at  /o,  has  been  suggested  to  improve  the  noise  figure  (NF)  of  an  FM  fink.  However, 
there  is  a  tradeoff  between  decreasing  the  dc,  which  decreases  shot  noise,  and  reducing 
the  signal  gain,  so  an  optimal  bias  point  must  be  found.  The  filter  cannot  be  biased 
exactly  at  the  null  or  the  fink  would  have  zero  output  current,  since  I  find  in  (2.5.28) 
that  the  output  is  proportional  to  the  bias.  This  is  consistent  with  experience  with 
carrier  suppression  on  IM-DD  finks. 

The  noise  contribution  from  optical  phase  noise  is 

S-pn  (2.6.1) 

71 

To  find  the  gain  and  noise  figure  of  the  fink,  the  frequency  modulation  depth  should 
be  written  in  terms  of  the  FM  laser  modulation  parameters.  It  is  proportional  to 
the  modulation  efficiency  rj,  in  units  of  Hz/A,  typically  of  the  order  of  a  few  hundred 
MHz  per  mA.  The  peak  input  current,  iin  produces  an  rms  input  power  when 
delivered  to  a  load  resistance  Rioad- 

=  iviin)  ^  =  27]"^ Pin/ Rioad  (2.6.2) 

Assuming  equal  input  and  output  loads,  the  power  gain  of  the  fink  is 

G  (2.6.3) 

The  noise  figure  of  the  fink  is  given  by  the  ratios  of  the  input  and  output  SNRs, 
assuming  a  thermal  noise  limited  input.  I  assume  that  the  dominant  sources  of  noise 
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are  phase,  shot  and  thermal  noises  at  temperature  T^.  The  relative  intensity  noise  is 
suppressed  through  the  balanced  detection. 


NF 


=1  + 
=1  + 


1  ‘^Q'idcRload  +  ^RloaFdc'^  ^  ^ 

G  ^  CksTK 

1  Tq  +  2idcA‘^^ 

G  +  2tdcV^A^kBTl 


(2.6.4) 

(2.6.5) 


A  useful  question  is  whether  it  makes  sense  to  low  bias  the  filter  in  an  attempt 
to  improve  the  noise  figure.  The  answer  depends  on  whether  the  designer  is  limited 
by  optical  power  available  or  by  the  maximum  photocurrent  the  photodetectors  can 
handle.  In  the  shot  noise  limited  regime,  the  noise  figure  is  improved  by  low  biasing 
the  filter,  as  long  as  the  optical  power  is  increased  to  maintain  a  fixed  dc  photocurrent 


NFsn  ~1  + 


T 


R 


load 


Tq 


2idcr]‘^A‘^kBTK 


(2.6.6) 


For  a  fixed  current,  for  which  the  optical  power  is  increased  to  maintain,  the  derivative 
of  the  NF  with  respect  to  the  bias  is 


dNF _  1  q 

The  noise  figure  is  always  improved  by  reducing  the  bias.  However,  the  phase  noise 
will  begin  to  dominate  over  the  shot  noise  when  Tq  <  2idcA‘^^,  and  any  NF  im¬ 
provement  will  be  negligible.  For  example,  with  Au  =  1  MHz  and  A  =  1/50  GHz  , 
choosing  a  bias  point  ^/T  <  0.5  only  makes  sense  if  the  maximum  dc  current  is  less 
than  160  fdA. 

If  the  available  optical  power  is  fixed,  idc  =  ^PoptT,  then  the  derivative  of  the  NF 
with  respect  to  the  bias  is  always  negative.  It  only  serves  to  reduce  the  gain  of  the 
fink  and  the  NF  gets  worse  with  the  lower  bias.  For  high  optical  powers,  the  NF  is 
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Figure  2.6.1:  Phase  noise  limited  noise  figure  versus  linewidth  and  modulation  effi¬ 
ciency 


Phase  noise  limited  NF  (dB) 


phase  noise  limited: 


NFpr.  ~1  +  Rload  0,  ^  (2.6.7) 

rj^kBTKTT 

This  is  independent  of  the  filter  bias  and  the  slope  of  the  filter.  Because  random 
frequency  fluctuations  are  added  to  the  optical  carrier  at  the  same  time  as  it  is 
modulated,  the  maximally  achievable  SNR  is  set  at  the  laser,  and  cannot  be  improved 
by  the  rest  of  the  system.  This  formula  sets  a  fundamental  relationship  between  the 
maximally  achievable  noise  figure,  the  laser  linewidth  and  the  modulation  efficiency. 
For  a  given  noise  figure  and  modulation  efficiency,  the  maximum  laser  linewidth  is 

Az/  =  {NF^  -  1)  (2.6.8) 

^load 

This  fundamental  relationship  between  modulation  efficiency,  linewidth  and  noise 
figure  is  plotted  in  2.6.1  for  a  range  of  typical  values. 


34 


2.6.2  Residual  intensity  modulation 


Residual  intensity  modulation  sets  a  lower  limit  on  the  distortion  for  an  ideal  balanced- 
detection  FM-DD  link.  The  effect  of  residual  IM  can  be  obtained  from  (2.1.9).  It  is 
difficult  to  write  a  general  expression,  but  it  is  possible  to  expand  some  individual 
terms.  In  lowest  polynomial  order  of  the  modulation  depth,  the  currents  of  interest 
are 


^dc^^Popt\ho,of  (2.6.9) 

if^  ^^PoptRe  {  [/3i  (hi,o/io,o  -  (2.6.10) 

+  ^mi  {hiflhlQ  +  /ro,o/i-i,o)  exp  [j27r/it]  | 

*2/1  I  (^2, 0^0,0  ~  2/ii^o^-i,o  +  ^o,o^-2,o)  (2.6.11) 

+mi/3i  [h2fihQQ  —  ho, 0^-2, o) 

exp  [j47r/it]}  (2.6.12) 

*2/i-/2  ^^Poptlne  {  [I3ll32  (-h2,-ih*  0  +  h2,ohl  i  +  ‘2hi,-ih*_,  ^  (2.6.13) 

—  ho,-ihl2,o  +  ho,ohl2,i) 


+  mil3i/32eP  (^2,0^0  —  h2,_ihg  q  -|-  ho-ih*_2  Q  —  ho,ohl2,i) 
+  2^2/3ie  (^2, 0^0,1  +  h2,-ihQ g  —  2hiflh*_i^^ 

—2hi^-ih*_iQ  +  hofih*_2^i  +  ho-ih*_2Q) 

+  rn‘l/32eP'^  (^—hi-ih*_i  Q  + 

+  -mim2/3i  (h2,ohg  -|-  hi  _ih*_-^^  Q  —  i  —  ho,-ihl2,o) 

exp  [j27r(2/i  -  /2)t]} 
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For  the  linear-field  filter,  these  can  be  written  as 

tdc  ^^PoptT  (2.6.14) 

if^  ^idcRe  {  [6f^2AT~'^P  +  mie^]  exp  [j27r/i  (t  -  r)]}  (2.6.15) 

*2/i~*rfc-Re|  +  mi6f^T~^PAeR  (2.6.16) 

+^2 1  exp  [j47r/i  (t  -  r)]  | 

*2/i-/2  ^idcRe  !^mi6f^6f^^A^T~^eR  +  m2Sj^^A^T~^e~R  (2.6.17) 

+  ^mlSAT~^PAeP'^  +  ^mim25/iT“^/^yl 

exp[j27r(2/i-/2)(t-r)]}  (2.6.18) 


Since  the  intensity  modulation  is  residual,  the  frequency  modulation  will  be  much 
greater  than  the  intensity  modulation.  Without  balanced  detection,  the  dominant 
second-harmonic  term  is  quadratic  in  the  FM,  and  the  dominant  IMD3  terms  are 
quadratic  in  FM  and  linear  in  IM.  I  assume  the  phase  difference  between  the  FM 
and  IM  is  0  =  0.  This  assumption  is  not  necessarily  correct  for  an  external  phase 
modulator,  for  which  the  PM  and  IM  should  be  in  phase,  with  0  =  — 7r/2.  For  a 
single  ended  detector,  the  distortion  currents  are  approximately 

*2/1  {Single)  ^idcT~^ cos  [dvr/i  {t  -  r)]  (2.6.19) 

*2/i-/2  {Single)  ^idcT~^A^^  [nii5f^5f^2  +  ni25)^)  (2.6.20) 

cos  [27r  (2/i  -  02)  {t  -  r)] 

With  balanced  detection,  the  dominant  second-harmonic  terms  are  linear  in  the  IM, 
and  terms  even  in  A  cancel  in  the  balanced  detection.  The  dominant  IMD3  terms  are 
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quadratic  in  the  IM,  and  terms  even  in  A  also  cancel.  The  currents  are  approximately 


*2/1  {Balanced)  ^idcT  ^^‘^midf^Acos  [47r/i  (t  —  r)]  (2.6.21) 

*2/i-/2  {Balanced)  ^idcT~^^‘^-A  (2mim25/i  +  (2.6.22) 

8 

cos  pTT  (2/i  -  /2)  {t  -  r)] 

It  is  useful  to  normalize  the  IM  to  the  FM.  One  method  of  normalization  is  to  look 
at  the  optical  power  the  IM  and  FM  contribute  to  the  first-order  optical  sidebands  in 
the  small  signal  approximation.  The  optical  power  in  the  first  order  sidebands  from 
the  IM  is  Poptnil/16.  The  optical  power  in  the  first  order  sidebands  due  to  the  FM  is 
Poptdj^/4:fi.  I  define  a  relative  residual  IM,  F,  by 

[rr^JlQ  mi  fi  r]^  f  ^ 

where  r]  is  the  FM  modulation  efficiency  and  r]m  the  IM  modulation  efficiency.  A 
second  quantity  that  should  be  normalized  is  the  RF  photocurrents  due  to  the  FM 
and  IM  after  the  discriminator.  Using  (2.6.15),  the  ratio  of  the  amplitude  of  the  two 
currents  is 


\if(FM)\  2ASf,  To 
\if{IM)\  ,n,VT  r 

where  I  define  a  discriminator  gain  parameter 


(2.6.24) 


(2.6.25) 
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Using  the  above  normalizations,  the  values  for  RIM  limited  OIP2  and  OIP3  are 


01 P2  {Single)  =32Rioadilc 

(2.6.26) 

OIP3{Single) 

(2.6.27) 

OIP2  {Balanced)  =2Rioadidc 

(2.6.28) 

OIP3{Balanced)  =-Rioadidc  f  ^  j 

(2.6.29) 

The  output  intercept  points  are  related  to  powers  of  the  ratio  of  RF  photocurrents 
caused  by  the  FM  and  IM.  For  FM-DD  links  to  exhibit  superior  third-order  and 
second-order  distortion  performance,  the  FM  laser  must  be  optimized  for  low  residual 
IM. 
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Chapter  3 


Numerical  link  analysis 


In  order  to  achieve  linear  discrimination  of  optical  FM,  I  desire  to  closely  approximate 
the  ideal  linear-field  optical  filter.  As  reviewed  by  [39],  optical  filters  can  be  designed 
using  digital  filter  design  techniques,  by  specifying  the  coefficients  of  the  z-transform 
representation  of  the  filter.  The  problem  of  discriminator  design  reduces  to  one  of 
choosing  the  best  coefficients.  A  physically-realized  optical  filter  is  limited  in  its 
number  of  stages  so  it  will  only  closely  approximate  the  desired  linear  frequency 
response. 

Finite  impulse  response  (FIR)  filters,  with  all  zeros  and  no  poles  in  their  z- 
transform  representations,  may  work  well  as  FM  discriminators  because  symmetric 
FIRs  can  be  designed  to  have  exactly  linear  phase.  In  this  chapter,  I  present  two 
sets  of  FIR  coefficients  chosen  using  the  maximally  linear  and  minimax  relative  error 
criteria  [40]. 

The  performance  of  the  resulting  filters  in  links  are  analyzed  using  the  small 
signal  model,  full  signal  model  and  a  numerical  simulation.  I  compare  the  linearity 
of  optical  frequency  discriminators  based  on  the  MRE  and  maximally  linear  criteria 
and  find  that  the  both  sets  of  filters  surpass  the  Mach  Zehnder  interferometer  (MZI) 
in  performance,  with  the  maximally  linear  filter  the  better  of  the  two.  I  also  present 
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Figure  3.1.1:  Numerical  model  of  an  FM  photonic  link  with  two  discriminator  filters 
and  balanced  detection 


i(t)  or  V(t) 

the  tradeoffs  between  filter  order,  third-order  nonlinearity  and  small  signal  gain. 

3.1  Numerical  link  simulation 

The  signal  at  the  output  of  the  link  is  simulated  by  creating  a  time  domain  waveform, 
e(t)  oc  exp[2nfct  +  2nri  f  i{t)],  performing  a  fast  Fourier  transfer  (FFT),  weighting 
the  frequency  domain  waveform  by  a  given  filter  transfer  function,  performing  an 
inverse  FFT  and  squaring  the  time  domain  waveform  to  obtain  the  photocurrent.  The 
simulation  process  is  illustrated  in  figure  3.1.1.  I  use  a  200  GHz  sampling  frequency 
and  a  10,000  point  FFT  to  convert  the  time  domain  waveform  to  the  frequency 
domain.  The  code  for  the  simulation  is  included  in  Section  A.l.  The  simulation  code 
includes  the  effects  of  residual  intensity  modulation  as  well  as  imperfect  common¬ 
mode  rejection  from  the  balanced  photodetection. 

3.2  Filter  coefficients 

The  frequency  discriminator’s  transfer  function  is  similar  to  that  of  a  digital  differ¬ 
entiator  filter:  a  ramp  with  constant  slope.  Two  sets  of  design  criteria  for  digital 
differentiators  known  in  the  literature  are  the  minimax  relative  error  (MRE)  criteria 
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and  the  maximally  linear  criteria  (maxlin).  The  MRE  criteria  minimize  the  maxi¬ 
mum  relative  error  over  a  chosen  band  of  frequencies.  The  maximally  linear  criteria 
fix  a  number  of  derivatives  of  the  transfer  function  at  a  chosen  frequency,  guarantee¬ 
ing  high  accuracy  around  a  small  frequency  band.  If  this  band  is  comparable  to  the 
bandwidth  of  our  modulation,  overall  I  expect  high  linearity. 

I  obtain  coefficients  for  MRE  filters  by  numerically  optimizing  at  half-band,  over 
the  range  of  normalized  frequencies  0.457r  to  O.SStt  (27r  is  the  full  ESR).  The  co¬ 
efficients  obtained  for  the  MRE  filters  are  shown  in  Table  3.2.1  on  page  44.  The 
maximally  linear  filter  is  also  optimized  at  half-band,  at  the  point  =  7r/2.  The 
coefficients  of  second,  sixth,  tenth  and  fourteenth-order,  maxlin  filters  are  shown  in 
Table  3.2.2  on  page  44  along  with  the  corresponding  coefficients  for  a  Mach  Zehnder 
interferometer.  The  filter  coefficients  are  illustrated  in  Eigure  3.2.1  on  page  42.  The 
transfer  functions  of  the  filters  are  compared  with  respect  to  the  ideal  one  in  Eigure 
3.2.2  on  page  43.  Eor  any  given  order,  the  maximally  linear  filter  appears  to  more 
closely  approximate  the  linear  discriminator.  This  example  does  not  preclude  the 
possibility  that  better  MRE  coefficients  could  be  found  using  a  different  numerical 
optimizations.  The  fact  that  the  MRE  filter  is  a  Type  III  linear  phase  EIR  filter 
(odd-length,  anti-symmetric)  and  the  maxlin  filter  is  a  Type  I  linear  phase  EIR  filter 
(odd-length,  symmetric)  also  affects  the  comparison. 

3.3  Comparision  of  models 

To  verify  the  analytical  small  signal  model  derived  in  this  work,  I  compare  the  small 
signal  model  (Section  A. 2)  with  the  large  signal  (Section  A. 3)  and  numerical  models 
(Section  A.l).  The  parameters  used  in  the  simulation  are  listed  in  Table  3.3.1  on  page 
45.  The  third-order  performance  of  a  link  with  an  MZI  is  shown  in  Eigure  3.3.1  on 
page  46.  The  three  models  are  consistent  up  through  modulation  powers  such  that 
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Figure  3.2.1:  Filter  coefficients  for  tenth-order  MRE  and  maximally  linear  filters 
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Figure  3.2.2:  Transfer  functions  for  MRE  and  maximally  linear  filters  compared  to 
ideal  discriminator  transfer  function 


Table  3.2.1:  Filter  coefficients  chosen  using  the  MRE  criteria  for  6th,  10th  and  14th 
order  filters  in  z-transform  representation 


Order 

6 

10 

14 

aO 

0.0318 

0.0060306 

0.0012721 

al 

-0.0806 

-0.01359 

-0.0027526 

a2 

0.2818 

0.0493 

0.012216 

a3 

0 

-0.1068 

-0.02426 

a4 

-0.2818 

0.2933 

0.060268 

a5 

0.0806 

0 

-0.1199 

a6 

-0.0318 

-0.2933 

0.2993 

a7 

0.1068 

0 

a8 

-0.0493 

-0.2993 

a9 

0.01359 

0.1199 

alO 

-0.0060306 

-0.060268 

all 

0.02426 

al2 

-0.012216 

al3 

0.0027526 

al4 

-0.0012721 

Table  3.2.2:  Filter  coefficients  for  the  MZI  and  2nd,  6th,  10th  and  14th  order  maxlin 
filters  in  z-transform  representation 


Order 

MZI 

2 

6 

10 

14 

aO 

-1 

2 

-1 

2tt 

-i/'i 

(27r)23 

-3/5 

(27r)27 

-10/7 

(27r)2li 

al 

1 

2 

1 

2 

0 

0 

0 

a2 

-1 

-9 

-25/3 

-98/5 

2tt 

(27r)23 

(27r)2'^ 

(27r)2li 

a3 

1 

2 

0 

0 

a4 

-9 

-150 

-490/3 

(27r)23 

(27r)2l 

(27r)2li 

a5 

0 

1 

2 

0 

a6 

-1/3 

(27r)23 

-150 

(27r)27 

-2450 

(27r)2li 

a7 

0 

1 

2 

a8 

-25/3 

(27r)27 

-2450 

(27r)2li 

a9 

0 

0 

alO 

=3/5 

-490/3 

(27r)27 

(27r)2li 

all 

0 

al2 

-98/5 

(27r)2li 

al3 

0 

al4 

-10/7 

(27r)2li 
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Table  3.3.1:  Simulation  parameters 


Parameter 

Value 

fi  (GHz) 

1.9 

/2  (GHz) 

2.0 

Impedances  (ohms) 

50 

Modulation  efficiency  (GHz/V) 

200 

Laser  linewidth  (MHz) 

0.5 

Filter  free  spectral  range  (GHz) 

100 

Optical  power  before  filters  (mW) 

100 

Optical  power  incident  upon  each  detector  (mW) 

12.5 

Photodetector  responsivity  (A/W) 

0.8 

Total  DG  photocurrent  (mA) 

20 

Third-order  IMD  (GHz) 

1.8 

Second  harmonic  (GHz) 

3.8 

the  fundamental  output  starts  to  saturate. 

The  third-order  performance  of  a  link  with  a  tenth-order  maxlin  filter  is  shown 
in  Figure  3.3.2  on  page  46.  The  three  models  are  consistent  for  small-modulation 
depths.  For  small-modulation  depths,  the  maxlin  filter  has  IMD3  several  orders  of 
magnitude  less  than  the  MZL  However,  the  small-signal  model  does  not  capture  the 
fast  increase  in  power  measured  at  frequency  2/i  —  /2  at  larger  modulation  depths. 
The  effect  is  due  to  increasing  amounts  of  optical  power  in  the  sidebands  that  fall 
outside  of  the  portion  of  the  filter  that  is  optimized  for  linearity.  The  bandwidth 
occupied  by  the  modulation  is  approximately  given  by  Carson’s  rule.  The  MZI  fink 
does  not  produce  second-order  distortion.  With  balanced  detection,  the  maxlin  fink 
has  very  little  second-order  distortion.  Assuming  a  common-mode  rejection  ratio  of 
30  dB,  the  distortion  is  plotted  in  Figure  3.3.3  on  page  47. 


3.4  Comparison  of  filters 


Using  the  numerical  and  small  signal  models,  I  compare  the  distortion  performance 
for  the  two  different  filter  optimizations.  The  fundamental  and  third-order  distortion 
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Figure  3.3.1:  Third-order  performance  of  an  FM-DD  link  with  an  MZI  discriminator 

Mach-Zehnder  discriminated  link 


Input  RF  Power  (dBm) 


Figure  3.3.2:  Third-order  performance  of  an  FM-DD  link  with  a  maxlin  discriminator 


lOth-order  maxlin  discriminated  link 


Input  RF  Power  (dBm) 
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Figure  3.3.3:  Second-order  performance  of  an  FM-DD  link  with  a  maxlin  discrimina¬ 
tor 

lOth-order  maxlin  discriminated  link 


powers  are  plotted  in  Figure  3.4.1  on  page  48  and  the  second-order  distortion  in  Figure 
3.4.2  on  page  48.  The  noise  floor  is  calculated  using  the  small  signal  model,  assuming 
a  500  kHz  laser  linewidth.  In  general,  the  higher-order  Alters  have  less  distortion  than 
the  lower-order  Alters. 

The  distortion  of  links  using  either  type  of  Alter  increases  much  faster  than  30 
dB  /  decade  for  large  modulation  power.  This  can  be  explained  by  observing  that 
in  frequency  modulation,  the  frequency  deviation  of  the  carrier  increases  with  mod¬ 
ulation  depth,  so  more  optical  power  is  spread  into  higher  order  sidebands.  For  high 
modulation  depths,  most  of  the  optical  power  lies  outside  of  the  range  of  frequencies 
for  which  the  Alter  is  optimized,  creating  more  distortion  than  for  low  modulation 
depths. 

The  FSR  of  the  Alter  is  a  design  parameter,  which  is  determined  by  the  delay 
elements  in  the  implementation.  A  smaller  FSR  gives  a  larger  slope,  meaning  better 
conversion  efAciency  from  FM  to  AM  and  more  link  gain.  It  can  also  affect  the 
linearity.  Figure  3.4.3  on  page  49  shows  OIP3  from  the  small  signal  model  versus  the 
FSR  of  the  Alter  and  Figure  3.4.4  on  page  50  shows  the  tradeoff  between  OIP3  and 
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Figure  3.4.1:  Third-order  distortion  versus  filter  length  and  type 


Figure  3.4.2:  Second-order  distortion  versus  filter  length  and  type 
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Figure  3.4.3:  Third-order  distortion  versus  filter  free  spectral  range 


gain  for  the  maxlin  filters.  By  increasing  the  FSR,  the  third-order  distortion  decreases 
much  faster  than  the  gain,  thus  the  OIP3  improves.  Taking  advantage  of  this  effect 
is  limited  because  the  noise  figure  degradation  will  reduce  the  SFDR,  and  residual 
IM  affects  the  linearity  more  for  less  efficient  discrimination.  Figure  3.4.5  on  page  50 
shows  the  SFDR  versus  FSR  calculated  using  the  numerical  model  to  find  where  the 
noise  floor  and  IMD3  intercept.  Saturation  of  the  achievable  SFDR  is  clearly  evident. 
The  choice  of  filter  FSR  must  carefully  balance  these  design  issues. 

For  maximally  linear  filters  with  larger  FSRs,  the  third-order  distortion  decreases 
faster  than  the  gain  decreases,  so  OIP3  monotonically  increases.  Higher  order  filters 
have  increasingly  better  distortion  characteristics:  by  increasing  the  filter  order  by 
four,  the  OIP3  increases  by  20  dB  or  more.  The  maximally  linear  filters  are  much 
better  than  the  MRE  filters  for  a  given  filter  order.  The  14th-order  MRE  has  higher 
OIP3  than  the  sixth-order  maximally  linear  filter,  but  lower  OIP3  than  the  tenth- 
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Figure  3.4.4:  Tradeoff  between  small  signal  gain  and  0IP3 


Figure  3.4.5:  Third-order  spurious-free  dynamic-range  versus  filter  free  spectral  range 
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order  maximally  linear  filter.  Our  results  suggest  that  a  microwave-photonic  fink 
using  a  250  GHz  FSR,  tenth-order  maximally  linear  filter  as  a  frequency  discriminator 
provides  20  dB  improvement  in  SFDR  over  a  MZI  discriminator. 
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Chapter  4 


Filter  implementation 

4.1  Planar  lightwave  circuit  filters 

An  FM  discriminator  approximating  the  ideal  linear  field  response  can  be  constructed 
using  planar  lightwave  circuites  (PLC).  The  transform  function  of  an  FIR  filter  can 
be  realized  in  PLC  with  just  MZIs  and  directional  couplers.  One  implementation  of  a 
multi-stage  optical  FIR  filter  in  PLC  is  the  lattice  filter.  The  lattice  filter  architecture 
has  a  low-loss  passband  and  requires  only  N+1  couplers  for  an  Nth  order  filter,  which 
are  advantages  over  other  optical  filter  architectures.  The  lattice  filter  architecture  is 
shown  in  Figure  4.1.1  on  page  53,  indicating  for  each  stage  the  coupling  coefficients, 
designated  by  k,  and  the  phase  shifts,  designated  by  (p.  Each  stage  has  a  unit  delay, 
z~^ .  The  dashed  lines  indicate  additional  filter  stages  ommitted  from  the  figure. 
A  recursion  relation  exists  that  transforms  between  given  filter  coefficients  and  the 
corresponding  coupling  ratios  and  phase  shifts  [39].  The  recursion  relation  for  the 
lattice  filter  design  gives  1024  solutions  for  a  tenth-order  filter.  One  particular  solution 
is  fisted  in  Table  4.1.1  on  page  53. 

Up  to  tenth-order  FIR  lattice  filters  have  been  implemented  in  PLC  for  various 
applications.  A  research  group  at  NTT  laboratories  has  extensively  explored  tunable 
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Figure  4.1.1:  FIR  lattice  filter  architecture 


(pi  (p2  (plO 


Table  4.1.1:  Filter  phase  and  coupler  parameters  for  a  tenth-order  maxlin  discrimi¬ 
nator _ 


Phase  shift 

Value 

Coupling  ratio 

Value 

Tunable  coupler  phase 

0 

Kq 

0.674271 

0.607391 

^2 

—71 

Ki 

0.635460 

0.648224 

0 

0.837472 

0.414954 

994 

—71 

0.514751 

0.770645 

0 

Ki 

0.918512 

0.289487 

^6 

0 

K5 

0.538578 

0.746782 

^7 

—77 

Kg 

0.918384 

0.289721 

0 

«7 

0.515217 

0.770179 

^9 

—71 

Ks 

0.837319 

0.415161 

^10 

—71 

Kg 

0.635828 

0.647842 

^10 

0.325694 

0.963443 
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optical  FIR  lattice  filters.  Tunable  coupling  ratios  are  implemented  by  using  sym¬ 
metric  Mach-Zehnder  interferometers  with  thermal-optic  phase  shifters.  A  diagram 
of  a  tunable  FIR  filter  is  shown  in  Figure  4.1.2  on  page  55.  The  intended  application 
is  dispersion  compensation,  but  because  the  filters  are  tunable,  they  can  be  used  for 
any  filter  transfer  function  desired,  including  FM  discriminators. 

The  group  has  fabricated  eight-order  filters  in  silica  with  chromium  heaters,  with 
50  GHz  FSR  [41]  and  200  GHz  FSR  [42],  and  arrays  of  fifth-order  filters  with  50 
GHz  FSR  [43,  44,  45].  They  claim  control  of  the  individual  phase  shifters  to  accuracy 
better  than  O.Olvr  radians.  For  the  fifth-order  filters,  to  reduce  the  required  bias 
power  on  the  heaters,  they  use  a  phase-trimming  technique  that  involves  introducing 
heating  induced  stress.  They  have  also  proposed  a  100  GHz  FSR  filter  in  a  refiection 
configuration  to  double  its  effective  length  [46]. 

A  collaboration  between  Siemens,  University  of  Kiel  and  IBM  Research  Zurich  has 
implemented  the  same  architecture  on  a  more  compact  silicon  oxynitride  platform  [47]. 
The  applications  include  both  EDFA  gain  equalization  and  dispersion  compensation. 
They  have  demonstrated  sixth,  seventh  and  tenth-order  filters  with  100  GHz  FSR. 
The  collaboration  has  explored  a  number  of  adaptive  feedback  approaches  for  setting 
the  filter’s  phase  shifters  [48]. 

•  Optical  spectrum  analysis  [49]:  They  have  used  an  optical  spectrum  analyzer 
to  compare  the  amplified  spontaneous  emission  spectrum  to  a  desired  intensity 
profile.  A  computer  running  the  Levenberg-Marquart  optimization  algorithm  (a 
modified  Gauss-Newton  algorithm)  varies  the  power  to  the  phase  shifters  until 
the  desired  profile  is  obtained. 

•  Electrical  spectrum  monitoring  [50,  51,  52]:  ESM  is  another  feedback  approach, 
where  power  at  certain  frequencies  are  used  as  a  feedback  mechanism.  Pilot 
tones  or  knowledge  about  the  signaling  over  the  link  determines  the  optimal 
choice  of  electrical  filters. 
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Figure  4.1.2:  Tunable  PLC  FIR  lattice  filter  architecture 


•  Eye  opening  [53,  54,  55]:  An  adaptive  feedback  approach  for  digital  signals  looks 
at  an  eye  diagram  and  uses  a  Levenberg-Marquart  optimization  to  maximize 
the  eye  opening. 

•  LMSE  /  minimize  ISI  [56,  57,  58,  50,  52]:  Another  method  for  digital  signals 
uses  minimization  between  the  decision  and  signal  as  a  feedback  signal.  These 
methods  are  not  suitable  for  analog  links. 

Eor  simplicity  and  cost,  setting  the  filters  coefficients  without  using  a  feedback  system 
is  desired.  The  IBM  collaboration  has  developed  a  calibration  procedure  to  produce 
a  table  look-up  for  tunable  coupler  and  phase  shifter  responses  versus  applied  tuning 
power  [49].  The  technique  uses  the  OS  A  approach  to  iteratively  tune  all  couplers  to 
zero  cross  coupling.  There  is  a  procedure  to  individually  characterize  each  tunable 
coupler  and  asymmetric  MZI  by  measuring  the  output  power  versus  tuning.  The 
filter  then  can  be  set  to  a  pre-calculated  inverse  system.  An  alternative  calibration 
approach  is  given  by  the  NTT  group  in  [59]  that  does  not  require  a  feedback  loop. 
The  approach  uses  incoherent  light  to  characterize  each  tunable  coupler,  and  coherent 
light  to  characterize  the  asymmetric  Mach  Zehnders. 

4.2  Monte  Carlo  simulation 

I  perform  a  Monte-Carlo  simulation  to  study  the  effect  of  imperfect  tuning  on  the 
resulting  link  distortions  caused  by  the  filter.  Eor  a  given  maximum  phase  deviation,  a 
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random  phase,  uniformly  distributed  between  plus  and  minus  the  maximum  deviation, 
is  added  to  the  path  length  of  each  symmetric  and  assymetric  MZI.  After  each  phase 
shift  is  varied,  the  new  filter  coefficients  are  calculated,  and  the  resulting  OIP3  and 
OIP2  are  calculated.  The  simulation  is  performed  for  100  trials.  Then,  the  simulation 
is  repeated,  but  the  symmetric  and  assymetric  MZI  phase  shifts  are  separated  out  to 
determine  if  the  system  is  more  sensitive  to  one  or  the  other. 

Figure  4.2.1  on  page  57  and  Figure  4.2.2  on  page  58  show  that  OIP3  and  OIP2 
degrade  significantly  if  the  phases  are  not  controlled  to  within  0.01  radians  of  their 
designed  values.  OIP2  is  more  sensitive  to  phase  errors  than  OIP3  because  the  phase 
difference  between  the  two  photodetectors  determines  the  cancellation  of  the  second 
harmonics.  The  approach  to  setting  the  filter  parameters,  as  discussed  in  Section  4.1, 
needs  to  be  able  to  meet  this  required  precision. 
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Figure  4.2.1:  Monte  Carlo  simulation  on  third-order  distortion 
Both  coupler  and  delay  line  errors 
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Figure  4.2.2:  Monte  Carlo  simulation  on  second-order  distortion 
Both  coupler  and  delay  line  errors 
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Chapter  5 


Conclusions 


In  this  report,  I  have  presented  frequency-modulated  microwave  photonic  links  with 
direct  detection  as  an  approach  to  achieving  links  with  low  noise-figure  and  high  lin¬ 
earity.  In  Chapter  1,  I  have  comprehensively  reviewed  previous  work  on  FM  discrim¬ 
inators,  including  early  work  with  Mach  Zehnder  interferometers,  and  later  attempts 
to  linearize  the  finks  using  “linear”  FM  to  IM  discriminators.  I  pointed  out  that 
this  approach  is  Hawed  since  it  ignores  coherence  and  combines  an  instantaneous  fre¬ 
quency  view  of  the  modulation  with  a  Fourier  frequency  spectrum  view  of  the  filter. 
As  an  alternative,  I  proposed  balanced  FM  discriminator  filters  that  are  linear  in  field 
transmission  versus  frequency  to  obtain  highly  linear  links. 

In  Chapter  2, 1  have  derived  large-signal  and  small-signal  analytical  expressions  for 
the  response  of  a  filtered  FM  fink.  Using  the  small-signal  model,  I  derived  figures  of 
merit  for  the  signal  to  noise  ratio,  third-order  intercept  point,  second-order  intercept 
point,  and  spurious  free  dynamic  range  for  a  fink  with  any  arbitrary  filter.  The  theory 
has  been  shown  to  be  consistent  with  previous  work  on  MZIs.  A  filter  with  a  linear 
FM  to  IM  ramp  has  been  shown  to  generate  large  amounts  of  distortion,  while  a  filter 
with  a  linear  FM  to  AM  ramp  theoretically  would  not  produce  any  third-order  or 
higher-order  distortion. 
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For  a  link  with  linear  FM  to  AM  conversion,  expressions  for  the  gain  and  noise 
figure  have  been  derived.  Low  biasing  the  filter  is  shown  to  improve  the  noise  figure 
only  if  the  photodetector  is  limited  to  very  small  amounts  of  optical  power.  Simplified 
theory  on  the  effect  of  residual  IM  has  been  developed,  showing  that  with  balanced 
detection,  OIP3  and  OIP2  quadratically  depend  on  the  ratio  of  RF  photocurrent  due 
to  FM  to  the  RF  photocurrent  due  to  IM. 

In  Chapter  3,  I  have  presented  a  numerical  model  for  simulating  the  distortion 
of  a  filtered  FM  link.  The  numerical  model  was  shown  to  be  consistent  with  both 
the  small-signal  and  large-signal  models  derived  in  this  report.  The  numerical  model 
was  used  to  compare  two  FIR  filter  approximations  to  the  linear  FM  to  AM  filter. 
Of  the  minimax  relative  error  criteria  and  the  maximally  linear  criteria  for  choosing 
the  filter  coefficients,  the  maximally  linear  criteria  has  been  shown  to  produce  lower 
distortion  for  any  filter  order,  and  to  improve  upon  on  the  MZI  in  OIP3  and  SFDR 
by  several  orders  of  magnitude. 

Finally,  in  Chapter  4,  I  reviewed  implementations  of  optical  FIR  filters  using 
the  lattice  filter  architecture  in  planar  lightwave  circuits.  Up  to  tenth-order  tunable 
arbitrary  filters  have  been  constructed  on  silica  and  silcon  oxynitride  platforms.  I 
have  reviewed  approaches  to  setting  the  filters’  phase  parameters  presented  in  the 
literature.  I  have  designed  an  implementation  of  a  tenth-order  maxlin  FIR  filter  us¬ 
ing  an  optical  lattice  filter.  I  simulated  its  sensitivity  to  the  phase  variations  on  its 
component  switches  and  delay  lines  using  a  Monte  Carlo  simulation  and  found  that 
the  phases  need  to  be  controlled  to  better  than  0.01  radians.  This  analysis  provides 
guideance  for  implementing  a  low-distortion  microwave-photonic  link  using  frequency 
modulation  and  discriminator  demodulation  with  a  tunable  FIR  filter.  Further  re¬ 
search  will  implement  the  filter  using  a  silica  PLC  device  and  measure  its  distortion 
performance  in  an  experimental  link. 


60 


Appendix  A 


Simulation  code 


A.l  Numerical  simulation  of  an  FM  link 

The  numerical  model  for  an  FM  link  is  based  on  the  model  in  Figure  3.1.1  on  page 
40.  It  includes  the  effects  of  residual  intensity  modulation  and  includes  the  effect  of 
the  common-mode  rejection  ratio. 

function  output  =  model_fft  (param ,  HaF,  HbF,  K,  fcMHz ,  PindBm,  CMRR) 
%MODEL_FFT  Numerical  model  for  a  frequency  modulated  link 
%  Input  param  =  [fsMHz,  N,  etaMHz,  G,  PmW,  R,  Rohm]; 

%  Input  Ha,  Hb  are  the  filter  transfer  functions  ,  over  the  range  of 
%  frequency  offsets  f  =  (—N/2:N/2—l)*fsMHz/N; 

%  Input  K  is  the  optical  power  split  [Pa/P  Pb/P]  between  the  two  filters 

%  Input  fcMHz  are  the  modulation  frequencies  in  MHz  [fl,  f2  ,  ■■■] 

%  Input  PindBm  are  the  modulation  powers  in  dBm  [PI,  P2 ,  ■■■] 

%  Input  CMRR  is  the  common  mode  rejection  ratio  for  balanced  detection 

%  in  positive  decibels  .  Matlab  accepts  ’Inf’  for  perfect  CMRR. 

%  Output  is  power  [PI,  P2 ,  ■■■]  in  dBm  at  frequencies 

%  [fl  ,  f2  ,  fn,  2*fl,  2*f2,  ...  2*fn,  2*fl-f2  ,  2*f2-fl] 

%%  Modulation  %% 
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fsMHz 

N 

etaMHz 

G 

RnW 

R 

Rohm 


=  param  ( 1 )  ;  %  FFT  sampling  frequency 
=  param  (  2  )  ;  %  FFT  number  of  points 

=  param  (3);  %  Frequency  modulation  efficiency  in  MHz/mA 
=  param(4);  %  RIM  normalization  parameter 
=  param  ( 5 )  ;  %  Optical  power  before  the  filters 
=  param(6);  %  Responsivity  of  photodetectors 
=  param  ( 7 )  ;  %  Input  and  output  load  impedance 


%  Monitor  the  output  for  second  and  third  order  distortions 

monMHz  =  [fcMHz*2,  abs  (fcMHz  (1)  *2 -fcMHz  ( 2 ) )  ,  abs  (fcMHz  (2)  *2 -fcMHz  ( 1 ) )  ]  ; 


%  C alculations  for  modulation 
PinW  =  1 0 .  "  ( ( PindBm  — 30)  / 10) ; 
lA  =  sqrt  (2*PmW/Rohm) ; 

ImA  =  IA*1000; 

df  =  etaMHz  *ImA ; 

b  =  df./fcMHz; 

m  =  G*b  *  2 ; 


%  Input  RMS  power  at  each  tone  in  Watts 
%  Input  peak  current  in  Amps 
%  Peak  amplitudes  for  each  tone  in  mA 
%  Frequency  modulation  depths 
%  Phase  modulation  depths 
%  Intensity  modulation  depths 


%  Generate  the  time  domain  signal 

linT  =  0;  XinT  =  0;  YinT  =  0;  %  Current,  FM,  and  IM 
for  X  =  1 :  length  (fcMHz) 

linT  =  linT  +  ImA(x)  / 1000*  sin  ( 2  *  pi  *fcMHz  (x )  *  ( 0  :N— l)/fsMHz  ) ; 
XinT  =  XinT  +  b  (x)  *  cos  ( 2  *  pi*fcMHz  (x )  *  ( 0  :N— l)/fsMHz  ) ; 

YinT  =  YinT  +  m(x)  *  sin  ( 2  *  pi*fcMHz  (x  )  *  (  0  :N— l)/fsMHz  )  ; 

end 


%  Generate  the  modulated  optical  signal  centered  at  fc 
EmodT  =  sqrt  (2  *PmW/ 1000)  *  sqrt  (1+YinT  ) .  *  exp  ( j  *XinT  )  ; 

%%  Filtering  %% 

EmodF  =  fft  shift  (  fft  (EmodT, N) ) /N;  %  Convert  to  the  frequency  domain 
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EoutaF  =  EmodF.  *HaF*sqrt  (K(  1 ) )  ;  %  Split  and  Filter 

EoutbF  =  EmodF.  *HbF*sqrt  (K(  2  ) )  ;  %  Split  and  Filter 

EoutaT  =  i  f  ft  (  i  f  f  t  s  h  i  f  t  ( EoutaF  )  ,N)  *N;  %  Convert  to  time  domain 

EoutbT  =  i  f  ft  (  i  f  f  t  s  h  i  f  t  ( EoutbF  )  ,N)  *N;  %  Convert  to  time  domain 

%%  D eteetion  %% 

%  Convert  baek  to  eleetrieal  domain 
loutaT  =  R*EoutaT  .*  conj  ( EoutaT  )/ 2  ; 
loutbT  =  R*EoutbT  .*  conj  (EoutbT  )/ 2  ; 

loutT  =  loutaT  —  loutbT  +  (loutaT  +  loutbT )  *10"  ( — CMRR.  /  1  0 ) ; 

%  Convert  the  eurrents  to  the  frequeney  domain 
liuF  =  fftshift  (  fft  (liuT  ,N))/N; 
loutF  =  fftshift  (  fft  (loutT  ,N))/N; 

FElec  =  (0:N/2-l)*fsMHz/N; 

%  This  gives  the  RMS  power:  multiply  by  2  for  two  sided 
%  speetrum  ,  square  it,  then  divide  by  2  for  RMS 
%  The  DC  power  should  not  be  multiplied  by  2. 

PiuF  =  2*Rohm*abs  ( liuF  (N/2+ 1  :end  ) )  .  "  2  ; 

PoutF  =  2*Rohm*abs  ( loutF  (N/2+ 1  :end  ) )  .  "  2  ; 

PinF(l)  =  PinF(l)/2;  PoutF  (1)  =  PoutF  (l)/2; 

PiuFdBm  =  10* loglO  ( PiuF )  +  30; 

PoutFdBm  =  10* loglO  ( PoutF )  +  30; 

%%  Analysis  %% 

%  Caleulate  the  output  power  for  eaeh  original  tone  in  dBm 
PoutdBm  =  zeros  (1,  length  (fcMHz  ))  ; 
for  X  =  1 :  length  (fcMHz) 

PoutdBm (x)  =  PoutFdBm (  find  (  FElec  fcMHz  ( x  )))  ; 
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end 


%  Calculate  the  output  power  at  other  specified  frequencies  in  dBm 
PmondBm  =  zeros  ( 1  ,  length  (monMHz)  )  ; 
for  X  =  1 :  length  (monMHz) 

PmondBm(x)  =  PoutFdBm(  find  (  FElec  monMHz ( x  )))  ; 

end 

%  Function  output 

output  =  [PoutdBm,  PmondBm]; 


A. 2  Small  signal  simulation  of  an  FM  link 

The  small  signal  model  for  an  FM  link  is  based  on  the  expression  in  (2.1.10).  It  does 
not  include  the  effects  of  residual  intensity  modulation,  but  includes  the  effect  of  the 
common-mode  rejection  ratio. 

function  output  =  model_smsig  ( param ,  HaF,  HbF,  K,  fcMHz ,  PindBm,  CMRR) 
%MODEL_SMSIG  Small  signal  model  for  a  frequency  modulated  link 
%  Input  param  =  [fsMHz,  N,  etaMHz,  G,  PmW,  R,  Rohm]; 

%  Input  Ha,  Hb  are  the  filter  transfer  functions  ,  over  the  range  of 
%  frequency  offsets  f  =  (—N/2:N/2—l)*fsMHz/N; 

%  Input  K  is  the  optical  power  split  [Pa/P  Pb/P]  between  the  two  filters 

%  Input  fcMHz  are  the  modulation  frequencies  in  MHz  [fl,  f2  ,  ■■■] 

%  Input  PindBm  are  the  modulation  powers  in  dBm  [PI,  P2 ,  ■■■] 

%  Input  CMRR  is  the  common  mode  rejection  ratio  for  balanced  detection 

%  in  positive  decibels.  Matlab  accepts  ’Inf’  for  perfect  CMRR. 

%  Output  is  power  [PI,  P2 ,  ...]  in  dBm  at  frequencies 

%  [fl  ,  f2  ,  ....  fn,  2*fl,  2*f2,  ...  2*fn,  2*fl-f2  ,  2*f2-fl] 

%%  Modulation  %% 
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fsMHz 


=  param  ( 1 )  ;  %  FFT  sampling  frequency 
N  =  param  (  2  )  ;  %  FFT  number  of  points 

etaMHz  =  param  (3);  %  Frequency  modulation  efficiency  in  MHz/mA 

G  =  param(4);  %  RIM  normalization  parameter 

RnW  =  param  ( 5 )  ;  %  Optical  power  before  the  filters 

R  =  param(6);  %  Responsivity  of  photodetectors 

Rohm  =  param  ( 7 )  ;  %  Input  and  output  load  impedance 

%  C alculations  for  modulation 
PinW  =  1 0 .  "  ( ( PindBm  — 30)  / 10) ; 
lA  =  sqrt  (2*PmW/Rohm) ; 

ImA  =  IA*1000; 

df  =  etaMHz  *ImA ; 

b  =  df./fcMHz; 

m  =  G*b  *  2 ; 

PoptW  =PmW/1000; 

%  Normalize  frequencies 
NormF  =  vectorize  (  inline  ( ’f*N/fsMHz^+^N/2+l  ’  ,  ’f’,  ’fsMHz’,  ’N’)); 

%  Calculate  the  output  power  for  each  original  tone  in  dBm 
n  =  1;  Xla  =  0;  Xlb  =  0; 
for  g  =  0:n 

Xla  =  Xla  +  ( — 1) " g/ f ac t o r i al (n— g) / f a c t o r i al ( g ) . *  ... 

HaF  ( NormF  ( ( n-g )  *  fcMHz  ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HaF(NormF(-g*fcMHz  , fsMHz  ,N) ) )  ; 

Xlb  =  Xlb  +  ( — 1) " g/ f ac t o r i al (n— g) / f a c t o r i al ( g ) . *  ... 

HbF  ( NormF  ( ( n-g )  *  fcMHz  ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HbF(NormF(-g*fcMHz  , fsMHz  ,N) ) )  ; 

end 

XI  =  Xla*K(l)  -  Xlb*K(2)  +  (Xla*K(l)  +  Xlb*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 

Pout  =  2/ (2  "  ( 2 *n ) )  *Rohm*R"2*PoptW .  "  2  .  * b .  "  ( 2  * n ) .  *  abs  (XI ) .  "  2  ; 

PoutdBm  =  10* loglO  ( Pout )  +  30; 


%  Input  RMS  power  at  each  tone  in  Watts 
%  Input  peak  current  in  Amps 
%  Peak  amplitudes  for  each  tone  in  mA 
%  Frequency  modulation  depths 
%  Phase  modulation  depths 
%  Intensity  modulation  depths 
%  Optical  power  in  W 


65 


%  Calculate  the  output  power  for  the  second  harmonics  in  dBm 
n  =  2;  X2a  =  0;  X2b  =  0; 
for  g  =  0:n 

X2a  =  X2a  +  ( — 1) " g/ f ac t o r i al (n— g) / f a c t o r i al ( g ) . *  ... 

HaF  ( NormF  ( ( n-g )  *  fcMHz  ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HaF(NormF(-g*fcMHz  ,fsMHz  ,N) ) )  ; 

X2b  =  X2b  +  (  — 1) " g/ fact  or ial (n—g) / fact o ri al ( g ). *  ... 

HbF  ( NormF  ( ( n-g )  *  fcMHz  ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HbF(NormF(-g*fcMHz  , fsMHz  ,N) ) )  ; 

end 

X2  =  X2a*K(l)  -  X2b*K(2)  +  (X2a*K(l)  +  X2b*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 

Pout 2  =  2/ (2  "  ( 2 *n ) )  *Rohm*R"2*PoptW .  "  2  .  * b .  "  ( 2  * n ) .  *  abs  (X2 ) .  "  2  ; 

PoutdBm2  =  10* loglO  ( Pout2 )  +  30; 

%  Calculate  the  output  power  for  the  third  order  distortion  in  dBm 

n  =  2 ;  p  =  1 ;  X3a  =  0 ;  X3b  =  0 ; 

fcMHzRR  =  fcMHz  ( [ 2  1]);  dfRR=df([2  1]); 

for  g  =  0:n 

for  k  =  0:p 

X3a  =  X3a  +  (—!)"( p+g+k ) / factorial ( n— g )/factorial(g)  ... 

/ factorial ( p— k )/factorial(k).*  ... 

HaF  ( NormF  ( ( n-g )  *  fcMHz+  (-p+k )  *  fcMHzRR ,  fsMHz, N)).*  ... 

conj  (HaF(NormF(-g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

X3b  =  X3b  +  (—!)"( p+g+k )  /  factorial  ( n— g  )/factorial(g)  ... 

/ factorial ( p— k )/factorial(k).*  ... 

HbF  ( NormF  ( ( n-g )  *  fcMHz+  (-p+k )  *  fcMHzRR ,  fsMHz, N)).*  ... 

conj  (HbF(NormF(-g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

end 

end 

X3  =  X3a*K(l)  -  X3b*K(2)  +  (X3a*K(l)  +  X3b*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 

Pout 3  =  2/ ( 2  "  ( 2 *  (n+p ) ) )  *Rohm*R"2*PoptW  .  "  2  .  *  . . . 

b .  ^  ( 2  *  n  ) .  *  (  dfRR .  /fcMHzRR) .  ^  ( 2  *  p  ) .  *  abs  (X3 ) .  ^  2  ; 
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PoutdBm3  =  10* loglO  ( Pouts )  +  30; 


%  Function  output 

output  =  [PoutdBm,  PoutdBm2 ,  PoutdBmS]; 

A. 3  Large  signal  simulation  of  an  FM  link 

The  large  signal  model  for  an  FM  link  is  based  on  the  expression  in  (2.1.9).  It  does 
not  include  the  effects  of  residual  intensity  modulation,  but  includes  the  effect  of  the 
common-mode  rejection  ratio.  Chosing  to  terminate  the  infinite  series  after  a  large 
number  of  terms  gives  a  result  with  arbitrarily  good  precision. 

function  output  =  model_lgsig  (param ,  HaF,  HbF,  K,  fcMHz ,  PindBm,  CMRR) 
%MODEL_LGSIG  Large  signal  model  for  a  frequency  modulated  link 
%  Input  param  =  [fsMHz,  N,  etaMHz,  G,  PmW,  R,  Rohm]; 

%  Input  Ha,  Hb  are  the  filter  transfer  functions  ,  over  the  range  of 
%  frequency  offsets  f  =  (—N/ 2  :N/2— l)*fsMHz/N; 

%  Input  K  is  the  optical  power  split  [Pa/P  Pb/P]  between  the  two  filters 
%  Input  fcMHz  are  the  modulation  frequencies  in  MHz  [fl,  f2  ,  ■■■] 

%  Input  PindBm  are  the  modulation  powers  in  dBm  [PI,  P2 ,  ■■■] 

%  Input  GMRR  is  the  common  mode  rejection  ratio  for  balanced  detection 
%  in  positive  decibels  .  Matlab  accepts  ’Inf’  for  perfect  GMRR. 

%  Output  is  power  [PI,  P2 ,  ■■■]  in  dBm  at  frequencies 

%  [fl  ,  f2  ,  fn,  2*fl,  2*f2,  ...  2*fn,  2*fl-f2  ,  2*f2-fl] 

%%  Modulation  %% 

fsMHz  =  param  ( 1 )  ;  %  FFT  sampling  frequency 

N  =  param  (  2  )  ;  %  FFT  number  of  points 

etaMHz  =  param  (3);  %  Frequency  modulation  efficiency  in  MHz/mA 
G  =  param(4);  %  RIM  normalization  parameter 

RnW  =  param  ( 5 )  ;  %  Optical  power  before  the  filters 
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R  =  param(6);  %  Responsivity  of  photodetectors 

Rohm  =  param  ( 7 )  ;  %  Input  and  output  load  impedance 

%  C alculations  for  modulation 
PinW  =  1 0 .  "  ( ( PindBm  — 30)  / 10) ; 
lA  =  sqrt  (2*PmW/Rohm) ; 

ImA  =  IA*1000; 

df  =  etaMHz*ImA ; 

b  =  df./fcMHz; 

m  =  G*b  *  2 ; 

PoptW  =PmW/1000; 

%  Normalize  frequencies 
NormF  =  vectorize  (  inline  ( ’f*N/fsMHz^+^N/2+l  ’  ,  ’f’,  ’fsMHz’,  ’N’)); 

%  Calculate  the  output  power  for  each  original  tone  in  dBm 

n  =  1;  p  =  0;  Xla  =  0;  Xlb  =  0; 

fcMHzRR  =  fcMHz  ( [ 2  1]);  dfRR=df([2  1]); 

for  g  =  -10:11 

for  k  =  —10:11 

Xla  =  Xla  +  ... 

b  es  s  e  1  j  (n+g  ,  df  .  /  fcMHz  ).*besselj  (k,  dfRR .  /  fcMHzRR) .  *  ... 

besselj  (g,df./  fcMHz  ).*besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 

HaF  ( NormF  ( ( n+g )  *  fcMHz+k *fcMHzRR ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HaF(NormF(g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

Xlb  =  Xlb  +  .  .  . 

besselj  (n+g  ,  df  .  /  fcMHz  ).*  besselj  (k,  dfRR .  /  fcMHzRR) .  *  ... 

besselj  (g,df./  fcMHz  ).*  besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 

HbF  ( NormF  ( ( n+g )  *  fcMHz+k *fcMHzRR ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HbF(NormF(g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

end 

end 

XI  =  Xla*K(l)  -  Xlb*K(2)  +  (Xla*K(l)  +  Xlb*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 


%  Input  RMS  power  at  each  tone  in  Watts 
%  Input  peak  current  in  Amps 
%  Peak  amplitudes  for  each  tone  in  mA 
%  Frequency  modulation  depths 
%  Phase  modulation  depths 
%  Intensity  modulation  depths 
%  Optical  power  in  W 
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Pout 


=  2*Rohm*R"2*PoptW .  "  2  .  *  abs  (XI ) .  "  2  ; 

PoutdBm  =  10* loglO  ( Pout )  +  30; 

%  Calculate  the  output  power  for  the  second  harmonics  in  dBm 
n  =  2;  p  =  0;  X2a  =  0;  X2b  =  0; 
for  g  =  -10:11 

for  k  =  —10:11 

X2a  =  X2a  +  ... 

b  es  s  e  1  j  (n+g  ,  df  .  /  fcMHz  ).*besselj  (k,  dfRR .  /  fcMHzRR) .  *  ... 

besselj  (g,df./  fcMHz  ).*besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 

HaF  ( NormF  ( ( n+g )  *  fcMHz+k *fcMHzRR ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HaF(NormF(g*fcMHz+k*fcMHzRR,fsMHz  ,N) ) )  ; 

X2b  =  X2b  +  .  .  . 

besselj  (n+g  ,  df  .  /  fcMHz  ).*  besselj  (k,  dfRR .  /  fcMHzRR) .  *  ... 

besselj  (g,df./  fcMHz  ).*  besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 

HbF  ( NormF  ( ( n+g )  *  fcMHz+k *fcMHzRR ,  fsMHz  ,  N ) ) .  *  ... 

conj  (HbF(NormF(g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

end 

end 

X2  =  X2a*K(l)  -  X2b*K(2)  +  (X2a*K(l)  +  X2b*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 

Pout  2  =  2*Rohm*R"2*PoptW .  "  2  .  *  abs  (X2 ) .  "  2  ; 

PoutdBm2  =  10* loglO  ( Pout2 )  +  30; 

%  Calculate  the  output  power  for  the  third  order  distortion  in  dBm 
n  =  2 ;  p  =  1 ;  X3a  =  0 ;  X3b  =  0 ; 
for  g  =  -10:11 

for  k  =  —10:11 

X3a  =  X3a  +  ... 

besselj  (n+g  ,  df  .  / fcMHz  ) .  *  b e s  s e Ij  (— p+k  ,  dfRR .  / fcMHzRR) .  *  ... 

besselj  (g,df./  fcMHz  ).*  besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 

HaF  ( NormF  (( n+g  )*fcMHz+ (-p+k  )*fcMHzRR,  fsMHz,  N)).*  .  .  . 

conj  (HaF(NormF(g*fcMHz+k*fcMHzRR, fsMHz  ,N) ) )  ; 

X3b  =  X3b  +  .  .  . 
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b  es  s e  1  j  (n+g  ,  df  .  / fcMHz  ) .  *  b e s  s e Ij  (— p+k  ,  dfRR .  / fcMHzRR) .  * 
besselj  (g,df./  fcMHz  ).*besselj  (k,  dfRR .  /  fcMHzRR ) .  *  ... 


HbF(NormF((n+g)*fcMHz+(-p+k)*fcMHzRR,fsMHz,N)).*  .  .  . 
conj  (HbF(NormF(g*fcMHz+k*fcMHzRR,fsMHz  ,N) ) )  ; 

end 

end 

X3  =  X3a*K(l)  -  X3b*K(2)  +  (X3a*K(l)  +  X3b*K(2) )  *10  ^  ( -CMRR.  /  1  0 ) ; 

Pout 3  =  2*Rohm*R"2*PoptW .  "  2  .  *  abs  (X3) .  "  2  ; 

PoutdBm3  =  10* loglO  ( Pout3 )  +  30; 

%  Function  output 

output  =  [PoutdBm,  PoutdBm2 ,  PoutdBm3]; 
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